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(i) 


ABSTRACT 


Average liquid and gas phase mass transfer coefficients during 
bubble formation have been determined using the carbon dioxide-water 
and the Be eat rdven water Systems respectively. The value of the 
mass transfer coefficient for the CO2-water system was found to be 
about 47 cm./hr. at 25°C which is lower than the values reported for 
the bubble rise period. 

The gas phase mass transfer coefficient for an ammonia-nitrogen 
feed containing 19 percent by volume ammonia was found to increase with 
the bubble frequency. The range of frequencies studied was 3 - 10 
bubbles/second and the transfer coefficient ranged from 0.25 meactee! 
hr-cm’-atm, at 3 bubbles/sec. to 0.32 gm-moles/hr-cm*-atm. at 10 
bubbles/second. 

The geometry of the growing bubble was evaluated in terms of 
its surface area and volume as a function of time by taking high speed 
motion pictures of the bubble at 600 frames/second and subsequently 
using a numerical integration method to calculate the area and volume 
from the coordinates of a number of points on the bubble periphery. 

The results obtained were accurate to within + 2%. The surface area of 


the growing bubble could be approximated by an equation of the form, 
A=A+k't 
fo) 


where A and A, are the surface area of the bubble at time t and time 


zero respectively, and k' is an empirical constant. 
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ici) 
A theoretical model based on the penetration theory and which 

iaakes use of the experimentally observed area-time relationship gives a 
reasonable estimate of the mass transfer occurring under liquid phase 


controlling conditions during the bubble formation period. 
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CHAPTER 1 


INTRODUCTION 

Many industrial applications involve the use of gas-liquid 
contactors, such as distillation columns, absorption columns, and 
sparged reactors. A rational design of these units requires a knowl- 
edge of the rate of the transfer processes occurring, and to improve 
our ability to predict the transfer rate it is useful to have an under- 
standing of the fundamental phenomena involved. 

Gas-liquid contactors often involve dispersing a gas into a 
liquid in the form of bubbles. Transfer to or from the bubbles takes 
place in at least three recognized regions, namely bubble formation, 
bubble rise, and bubble breakup at the liquid surface. 

Although considerable experimental data and theoretical anal- 
yses are available for mass transfer during bubble rise, the analysis 
of transfer during formation has been rather limited. The only reported 
experimental results for mass transfer during bubble formation are 
those of Kalra [49 ] who studied the transfer from carbon dioxide bubbles 
growing in a moving water phase. His results indicate a higher rate of 
mass transfer Ck, = 100 to 225 cm/hr) occurring in this regime as com- 
pared to that obtained during bubble rise (k, = 125 cm/hr). In contrast, 
Calderbank and Patra [15] have reported no such enhanced rate based on 
a direct study of mass transfer for carbon dioxide, sulphur dioxide, 
and acetylene bubbles formed in a stagnant water phase. But they have 
not given any experimental results. Some other quantitative inform- 


ation has been obtained indirectly by Calderbank [9], Haselden et. al. 
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Estas Licht et. al. 576 58], and Dixon et. al. [27] by treating mass 
transfer during bubble formation as an end effect and has been inter- 
preted to indicate a significant amount (25 - 50%) of the total mass 
transfer occurring during bubble formation. These authors had extrap- 
olated the bubble rise data to zero liquid height and had assumed it 
to be applicable to bubble formation. 

The above mentioned studies, although limited, were for mass 
saree during bubble formation under liquid phase controlling transfer 
conditions. No investigation is believed to have been undertaken to 
date for mass transfer from a bubble during its growth when the resist- 
ance to transfer lies in the gas phase. However, Calderbank [9], 
Thorogood [77], and Haselden et. al. [37 ] have attempted to estimate 
the mass transfer into a growing bubble from a liquid phase under gas- 
phase controlled conditions and conclude that a considerable amount 
(20 - 50%) of the possible transfer may occur during formation of the 
bubble. 

The lack of published data and the presence of uncertainty as 
to the amount of mass transfer occurring in the bubble formation region 
justified the study undertaken here. Two systems, namely the carbon 
dioxide-water system and the ammonia-nitrogen-water system were invest- 
igated in a gas-liquid column by using an experimental technique similar 
in approach to the one described by Kalra [49]. The effect of varying 
the gas and the liquid flow rate on the transfer coefficient was 
studied. It is expected that the main resistance to mass transfer is, 
in the continuous phase for the C02 - water system, and in the dis- 
persed phase for the ammonia-nitrogen-water system. 


A review of the theoretical approaches to predicting mass 
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transfer during formation is given and in some cases modifications are 


proposed. 
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CHAPTER 2 


LITERATURE REVIEW 
To analyze mass transfer from gas bubbles during formation 
it is helpful to have information on the mechanics of bubble formation, 
and any experimental data available for mass transfer rates, and a 
theoretical understanding of the transfer process. Comprehensive 
reviews on bubble phenomenon have been undertaken by Jackson [46], 
Valentin [80], Calderbank [11], and Kalra [49]. 
In this literature review an attempt is made to summarize the 
conclusions of previous workers under the following categories: 
(1) Effect of variables on bubble behaviour 
(2) The bubble shape during growth 
(3) Mass transfer in the continuous phase 


(4) Mass transfer in the dispersed phase. 


2.1 Effect of Variables on Bubble Behaviour 
The principal factors which might be expected to affect the 

diameter of a bubble formed at an orifice are: 

(1) the volumetric flow-rate of gas through the orifice 

(2) the orifice diameter 

(3) the gas density and viscosity 

(4) the liquid density and viscosity 

(5) the surface tension 

(6) the wetting properties of the orifice material 


(7) the pressure drop across the orifice 
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(8) the chamber volume below the orifice 

(9) the submergence of the orifice 

(10) the shape of the orifice 

(11) the angle of inclination of the orifice; and 

(12) the hydrodynamics of the continuous phase. 

Some of these factors are clearly interdependent, for example, 

the orifice pressure drop is related to the gas flow, the orifice 
diameter, and the submergence. The effect of some of the above var- 


iables is discussed below. 


2.1.1 Effect of Gas Flow Rate 
As the gas flow rate is increased, four different bubble form- 
ation regimes are encountered, namely 
(a) constant volume region 
(b) slowly increasing volume region 
(c) constant frequency (laminar) region; and 
(d) turbulent region. 
The change from one region to the other has been characterized 


in terms of the orifice Reynolds number which is defined as, 


(a) Constant Volume Region 
At low gas flow rates (Re, upto 100 or 150), the gas bubble 
forms discretely and slowly at the orifice and is substantially con- 


stant in volume. Thus, the frequency of bubble formation is propor- 
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tional to the gas flow rate. The diameter of the bubble at detachment 
from the orifice can be calculated by a force balance between the 


surface tension and the buoyant forces as follows: 


Buoyancy force, 


F, = (1/6) d,* Apg = ------------------------------------------- (2.2) 


Surface tension force, 


F = ™D,y(COS 0.) | -------------------------------------------- (os } 


If the orifice is perfectly wetted by the liquid, then, 


Bueno) peya Vee ee a eee (2.4) 


Van Krevelin et. al. [s1], Quigley et. al. [69], Hughes et. al. 
[43], Coppock and Meikeljohn [19], Leibson et. al. [52] and Datta et. al. 


[22] are in agreement with the above observations. 


(b) Slowly Increasing Volume Region 

This is a transition region for the changeover from the con- 
stant volume to the constant frequency (laminar) region. A range of 
Re, of 100 - 300 has been suggested for this region by Valentin [80]. 
Neither the bubble volume nor the bubble frequency are constant in this 
region of bubble behaviour. Davidson and Shuler [25] have carried out 
experiments in this region and report an increasing predominance of the 
inertial effects. The following equation for bubble diameter is sug- 


gested, 


/4 


Man —-----------~--+----~-----~-- (2.5) 
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d, = const. (D¥/Ap) 
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10 
(c) Constant Frequency (Laminar) Region 


With a further increase in gas flow rate (Re, about 1,000 to 
2,100) the bubble frequency becomes independent of the gas flow rate 
and approaches a constant maximum value, which is a function of the 
orifice diameter. The bubble volume, however, increases with increas- 
ing gas rate. 

Davidson and Amick [24] have reported the value of this con- 
stant frequency as 15/sec. for an orifice diameter of 0.32 cm. and 
Eversole et. al. [28] found the value to be 45/sec. for a 0.02 cm. 
diameter orifice. The former authors found the following empirical 


relation to hold for the bubble diameter, 


bubble diameter, cm. 


z= 
an 
(i) 
ny 
@ 
Q. 
Ph 
i] 


e 
il 


gas flow rate to orifice, cm.°/sec. 


D_ = orifice diameter, cm. 


(d) Turbulent Region 

As the gas flow rate is increased further to a value of Re, 
greater than 2,100 a range of bubble sizes is formed, i.e. bubble 
break-up or coalescence starts occurring very near the orifice. Very 
small satellite bubbles of about 5 - 50 micron diameter are produced 
around the primary bubbles and cause considerable disturbance to the 
surface of the latter thus enhancing mass transfer. Leibson et. al. 


[52] and Calderbank [10] worked in the Reynold's number range of 
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2,100 - 10,000 and are in agreement with this behaviour. 
2.1.2 Effect of Orifice Characteristics 


(a) Orifice Diameter 

During the constant volume region the bubble size is propor- 
tional to the cube root of the orifice diameter as can be visualised 
from equation (2.4). The size of the bubble will eventually become 
smaller than the orifice diameter as the orifice size is increased. 
This causes unstable bubble formation and can result in liquid weeping 
through the orifice. Davidson and Amick [24] have found that orifices 
of radius larger than 0.7 cms. could not produce stable bubbles in the 
constant volume region. 

In the constant ae ene region Van Krevelin et. al. [81] have 
reported that the bubble volume is independent of orifice diameter. 
Davidson and Amick [24], however, found it to vary as 0.43 power of the 


orifice radius. 


(b) The Angle of Inclination of the Orifice 

Datta, Napier and Newitt [22] from their extensive experimental 
data found that the bubble size from an inclined orifice was smaller 
than from a vertical one. This could be explained by the tendency of 
the buoyancy forces to drag the bubble across the orifice, thus allow- 


ing the surface forces to operate over only a portion of it. 


(c) The Orifice Surface Wettability 


This factor is important for preventing liquid leakage through 
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the orifice. Stable bubble formation requires that the orifice surface 
be non-wettable by the liquid. Bowman et. al. [8] have used teflon 
covers on the orifice to improve the non-wettability for cases where 


liquid leakage occurred due to the orifice diameter being too large. 


2.1.3 Effect of Chamber Volume 

The importance of the chamber volume on bubble behaviour has 
been emphasized by Hughes et. al. [43], Davidson and Amick [24], Bowman 
and Johnson [8], Hayes et. al. [38], Leibson et. al. [52] and Mahoney 
et. al. [62]. The chamber volume is defined as the volume below the 
orifice to the nearest restriction in the gas line. The bubble sta- 
bility can be controlled effectively by varying the chamber volume 
which also affects the volume of bubble. A capacitance number has 


been defined by Hughes et. al. as, 


The formation of stable discrete bubbles occurs when No is 
below a certain critical value. For chamber volumes where No is 
greater than the critical value the bubble formation becomes unstable 
and a phenomena referred to as chain bubbling occurs. Davidson et. al. 
[24] and Hughes et. al. [43] found the bubble behaviour to be independ- 
ent of N. for No <0.8 at low gas flow rates and for No <0.2 at higher 
flow rates. It also becomes independent of No if very large chamber 
volumes are used (constant pressure condition). 


A different approach to the influence of chamber volume has 
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13 
been taken by Davidson and Shuler [25, 26], and more recently by Kumar 
and co-workers [50, 103 72]. two limiting cases of bubble formation 
ae defined wherein the chamber volume does not have any significant 
influence on the stability of bubble formation. They are; a) form- 
ation under constant flow conditions, and b) formation under constant 
pressure conditions. In the former case the flow rate of the gas into 
the bubble is maintained constant by causing a large pressure drop 
across the orifice, and in the latter the pressure of the gas below 
the orifice is maintained constant by using a chamber volume very much 
greater than the critical value defined in eqn. (2.7). However, it 
has been shown [25, 26, 50, 70, 72] that the dependence of bubble 
volume on other variables is different for the two cases. This may 
account for some of the discrepancies between the results of different 


workers on the effect of variables on bubble phenomena. 


2.1.4 Effect of Liquid Seal Height 

Khurana and Kumar [50] have recently investigated the effect 
of liquid seal height on bubble behaviour both theoretically and ex- 
perimentally for such chamber volumes that both the pressure inside 
the chamber and flow rate into the bubble are time dependent. The 
experiments were conducted for air bubbles in a water-alcohol mixture 
and air flow rates were varied from 3 cc/sec. to 20 cc/sec. through an 
orifice of diameter 0.3 cm. Orifice submergence was varied from 15.6 
cm. to 128 cm. and two chamber volumes of 65 and 600 cc were studied. 
They report that the bubble volume reduces exponentially with increase 


in orifice submergence. Bowman et. al. [8] and Leonard et. al. [54] 
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14 
have also shown a similar effect. 
However, Quigley et. al. [69], Davidson et. al. [24], Coppock 
and Meikeljohn [19], Hayes et. al. [38], and Towell et. al. [79] have 


all reported a negligible effect of seal height on bubble formation. 


2.1.5 Effect of Physical Properties of the Liquid 


Some workers have varied liquid properties considerably in 

their investigations, but their conclusions are largely contradictory. 

Datta et. al. [22] found that the bubble size decreased slowly 
with increasing liquid viscosity. On the other hand Quigley et. al. 
[69 ] observed a slow increase and Davidson et. al. [24] a marked in- 
crease. Coppock and Meikeljohn [19] found no effect of viscosity on 
bubble volume. 

An increase in liquid density was found to decrease the bubble 
size by Coppock et. al. [19], Benzing and Myers [4] and Davidson and 
Shuler [25]. Contrary to this Quigley et. al. [69] state that it has 
no effect. 

The bubble size was reported to increase with an increase in 
surface tension by Datta et. al. [22], Coppock and Meikeljohn [19] 
and Benzing and Myers [4]. However, Quigley et. al. [69] and Siemes 
et. al. [74] claim that it has no effect. This can perhaps be explain- 
ed from Davidson and Shuler's results [26] who found that surface 
tension has a negligible effect when the experiments are carried out 
under constant flow conditions and high pressure drop, but for constant 
pressure conditions the effect is appreciable. In other words, sur- 
face tension is hendeeeas only when inertial forces are relatively 


unimportant. 
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Delsow er fect of Liquid Velocity 
The experimental data of Sullivan et. al. [76] indicates that 
the liquid velocity past an orifice does not have any significant 
effect on the diameter of bubbles. Air bubbles at air-flow rates of 
0.5 - 100 cc/sec. through orifice diameters of 0.16 to 0.32 cm. mate | 
formed in water flowing at velocities of 0.34 to 2.5 cm./sec. horizon- 


tally to the orifice. 


.2 The Bubble Shape During Growth 
Kalra [49] employed high speed photography to determine the 


shape of the bubble during its growth from a horizontal orifice and 
his findings are depicted in Figure 1 which shows a typical bubble 
during its various stages of formation. The gas bubble oscillates in 
the axial direction during the initial stages of growth which is due 
to the break-off of the previous bubble. However, these oscillations 
cease fairly quickly and the bubble grows steadily and uniformly 
thereafter in an almost spherical pattern until the volume of the 
bubble approaches a limiting value. A neck then develops and the 
bubble begins oscillating in its axial direction, finally breaking off 
from the orifice. Leibson et. al. [52] have suggested this necking tu 
be due to inward liquid circulation at the level of the orifice. Cine 
pictures of a bubble growing from an orifice of 0.25 inch diameter as 
obtained by McCann and Price [60] follow the same general pattern as 
described above. 

The area of the growing bubble as a function of time is re- 


quired in order to determine the average mass transfer coefficient 
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during the entire formation regime and has been investigated by Kalra 
[49], and Calderbank and Patra [15]. Poutanen et. al. [68] and 
Halligan et. al. [34 | have developed equations to describe the bubble 
profile but these do not predict the variation with time. 

From the photographs of bubbles forming at an orifice Kalra 
[49] evaluated the surface area and volume at various Stages of growth 
by reading the coordinates of various points on the bubble periphery 
and subsequently using a computer oriented numerical integration ap- 
proach. His experimental results indicated the area of bubble upto 


neck formation to follow the relationship 


where k is a constant. 
Calderbank et. al. [15] through a similar frame-by-frame anal- 


ysis of the cine film records using an analytical projector report an 


empirical correlation of the type 


where t is the time from commencement of bubble formation, te is the 
total formation time, and AL is the maximum area obtained during the 
bubble growth. 


The volume of the bubble varies linearly with time provided 


the gas flow rate to the orifice is held constant [49]. 
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2.3 Mass Transfer in Continuous Phase 


2.3.1 Transfer During Bubble Rise 


2.3.1.1 Experimental Data: 

Mass transfer from a bubble rising singly in a liquid has been 
studied quite extensively and a considerable amount of experimental 
data are available which are summarized in Table 1 for the carbon 
dioxide-water system. Considerable scatter in the values of transfer 
coefficient exists and can be attributed partly to different experi- 
mental techniques used by various authors. The bubble rise velocity 
lies in the range of 20 - 40 cm./sec. depending on the bubble size. 

The value of the continuous phase coefficient reaches a max- 
imum as the bubble diameter increases but subsequently decreases with 
any further increase in bubble diameter. Valentin [80], Zieminski and 
Raymond [85], and Calderbank et. al. [16] report such a peak for equiv- 
alent bubble diameters in the range of 0.1 - 0.4 cm. which may be due 
to a progressive transition from rigid bubble behaviour to circulating 
bubble behaviour as the bubble diameter increases. The subsequent fall 
in the value of transfer coefficient with increase in diameter of the 
circulating bubble is due to an increase in the value of the contact 


time. 


2.3.1.2 Theoretical Analyses: 
Calderbank [11], Redfield and Houghton [71], and Valentin [80] 


have categorized the analysis of mass transfer during bubble rise as 
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20 
follows: 
(a) transfer from rigid spheres 
(b) transfer from circulating spheres. 

Small bubbles below about 0.1 cm. are considered to have no 
internal circulation and consequently do not have any velocities at the 
interface. On the other hand large bubbles have a mobile interface 
because of the fluid inside them being under constant circulation. 

This internal circulation is caused by viscous forces between the two 
fluid phases. 

Mass transfer from each of the above mentioned categories is 
generally analyzed for different magnitudes of the Reynolds number 


which is defined as, 


Ren= aL ST (2.10) 
L 
where de = equivalent diameter of bubble 
U = terminal velocity of rise of bubble 
PL = continuous phase density 
Wy = continuous phase viscosity. 


Solid Spheres: 

Levich [55] obtained the following equation for creeping flow 
around a solid sphere (Re <1) by assuming a thin concentration boundary 
layer (100 < Pe < 1000) in the continuous phase. 
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Friedlander [29] and Baird et. al. [2] have also obtained 


similar results with proportionality constants of 0.89 and 0.99 re- 
spectively instead of 1.01. 
The Stokes equation can be used for the bubble rise velocity 


in equation (2.11) and is, 


where V is the kinematic viscosity of the liquid in centipoise, V 
is the bubble volume in c.c., and U is in cm./sec. 

Under conditions where inertial effects can no longer be 
neglected (Re >>1) a semi-theoretical approach is usually used. 


Griffith [31] proposed the following relationship, 
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As Re approaches zero, the above equation reduces to 
Frosseling's equation for steady diffusion from a stagnant sphere which 


Les 


Circulating Spheres: 


For creeping flow (Re <1) around a circulating sphere when 
the assumption of a thin concentration boundary layer (Pe *?1) holds, 


Levich [55] obtained 
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A similar equation has been reported by Ward et. al. [84] for 
the case when the interface moves quickly with respect to the main 
stream velocity. However, when the interfacial velocity is a fraction 


of the main stream velocity, Ward et. al. find that 


Se Une e oe oe bewsen (555 doamd “Rene i) Sol eek he (2.16) 


The dividing line between equations (2.15) and (2.16) is at 
about Pe equal to 1,000. 
The Hadamard - Rybczynski equation for the terminal velocity 


in this regime is: 


U = 126 2/3 3 (2/17) 
Vv 
For cases where Re >>] Lochiel and Calderbank [59 ] assumed the 
thickness of the velocity boundary layer to be far less than the equiv- 
alent spherical radius of the dissolving bubble and the concentration 
boundary layer to be far thinner than the velocity boundary layer 
(Sc >>1). On this basis they derived: 
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The indicated range of Re for this equation is 100 to 400. 


When Re is very large equation (2.18) reduces to the Boussinesq 
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[7 ] equation for circulating spheres in potential flow, 


Sh = 1.13 Pee —--—----~---~+-~~~---~-~-~-------~-- +--+ +--+ ++ (2.19) 
an 
U 
ki = 1.13 Dy, 0.5 -------------------------------------------- (2.20) 
(>) 


The bubble rise velocity to be used in equations (2.19) and 


(2.20) is derived by Levich [55] and is, 


For gas bubbles in liquids Higbie [41] suggested that the mean 
exposure time of the elements of interface can be taken as the time 
for the bubbles to rise through a distance equal to their vertical 
diameter. With this assumption equation (2.20) becomes the familiar 


Higbie equation, 


2.3.2 Transfer During Bubble Formation 

The study of mass transfer during bubble formation has been 
rather limited. The only reported experimental results are those of 
Kalra [49] for carbon dioxide bubbles in water which indicate an en- 
hanced rate of mass teenerer (kK, = 100 to 225 cm./hr.) in comparison 
to that during bubble rise (k, = .125 cm./hr.) for the same equivalent 


bubble diameter. Calderbank and Patra [15] have reported that the 
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average liquid-phase mass transfer coefficient during the formation 

of a bubble at a submerged orifice is less than the value observed 
during its subsequent ascent but have not given any experimental re- 
sults. Some other quantitative information has been obtained indirect- 
ly by Calderbank [9], Haselden et. al. [37], icht fét real, [57, 58] and 
Dixon et. al. | Same} by treating mass transfer during bubble formation 

as an end effect and has been interpreted to indicate a significant 
amount (25 - 50%) of the total mass transfer occurring during bubble 
formation. These authors had extrapolated the bubble rise data to 

zero liquid height and had assumed it to be applicable to bubble form- 
ation. The only other experimental studies have been on mass transfer 
during drop formation. 

Popovich et. al. [67] have shown that the models proposed for 
mass transfer into the continuous phase from a growing drop by Licht 
and Pansing [58], Heertjes et. al. [40], Groothius and Kramers [32], 
Ilkovic [44, 45] and Michels [64] can all be expressed by the general 


equation 


N' = Bk, t eb (Aket1)/2 1 0p 9) CCP A oaas 


where the variation of drop area with time is given by 


and the variables are defined as: 


N' = mass transferred during bubble formation per unit area 
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ky = empirical constant 
ko = empirical constant 
te = formation time 

8 = empirical constant 


The value of 8 in equation (2.23) takes different values de- 
pending on the mechanism assumed. Popovich et. al. found their ex- 
perimental data on liquid drops to correspond to Ilkovic's model. The 
value of 8 is 1.31 for this model provided spherical drop growth is 

a 2. -2/3 
assumed. In other words when ky = Td te 


d is the drop diameter at the end of formation period. 


, and ko = 2/3, where 


However, the shape of bubbles forming at an orifice is generally 
less nearly spherical than liquid drops and hence the applicability of 
the analogy with liquid drops is uncertain. 

Calderbank and Patra [15] studied the mass transfer during bubble 
formation for carbon dioxide, acetylene and sulphur dioxide bubbles in 
water. They proposed an empirical correlation for the variation of 
bubble surface area with bine which was given as equation (2.9). It 
was assumed that each freshly formed element of surface is segregated 
from older elements. Using the theory for diffusion into a semi-infinite 
body Calderbank et. al. derived the following equation for the quantity 


of material dissolved during bubble formation: 


D 
L 
ed abivs6! (ci®> Ceyoas pe A Ge ah ae no roca can ee anar al TEES (2.25) 


The average mass transfer coefficient was calculated by assuming 


that the quantity of material transferred as given by equation (2.25) is 
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26 
the same as would be dissolved if the bubble was formed instantaneously 


and remained at its maximum size for a time ty. This implies 


From equations (2.25) and (2.26), 


All of the experimental results of Calderbank and Patra [15] 
with carbon dioxide and acetylene bubbles in distilled water agreed 
well with equation (2.27). 

Mass transfer coefficients during bubble formation and bubble 


rise can be compared from equations (2.20) and (2.27). One obtains 


d 
k, (formation) . 0.78/(-£-) site ak a 0 ici op po stra ta EE (2528) 
k (rise) Ute 


If 'g' is assumed to be the vertical distance between rising 


bubble centres, then for discrete bubble formation from the orifice, 


and equation (2.28) becomes, 
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Hence the mass transfer coefficient during bubble formation 
as defined by equation (2.26) is always lower than that during bubble 
rise as (d,/%) is less than one if there is no bubble coalescence at 
the orifice. 

Beek and Kramers [3] have analyzed the mechanism of mass 
transfer from an expanding spherical surface when the resistance to 
mass transfer lies in the continuous phase. They assumed that circula- 
tion existed inside the growing surface and that the increase in sur- 
face occurred by stretching of the surface. For the case where the 
depth of penetration could be considered small compared to the radius 
of the cavity, they derived the following expression for the amount 


transferred during formation: 


“1 See 6 C, eas set ASL ~--------- (2.31) 
f ] 
(7+4p) 1? 
Aes ee 
when _P >>1, and the radius of the bubble at any time is 
D3 
expressed as, 
R(t) = (3 cas as mr rrr (2.32) 


2.4 Mass Transfer in Dispersed Phase 


2.4.1 General Models 
The mechanism of internal mass transfer has been discussed by 
Newman [66 ], Kronig and Brink [51 ], Crank [20], Harriott [36], 


Skelland and Wellek [75], Miller [65], Johns et. al. [47], Thorogood 
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[77], and Calderbank [11]. 


The simplest model is Newman's [66] who studied the drying of 
porous solids, spherical in shape, during that portion of the drying 
time in which the surfaces are substantially dry and evaporation takes 
place as rapidly as liquid can diffuse from the interior of the sur- 
faces. He expressed the differential equation for unsteady state 
diffusion across a shell of thickness, dr, within the sphere at radial 


position, r, as 


a) 


Beep (We + 
or2 r 


ot G 


Kit 


where Da is the diffusivity in the dispersed phase, and c is the con- 
centration at radial position, r, at time, t. 
The following boundary conditions were specified: 
(a) the initial liquid concentration is uniform throughout the 
e0lid.51,6.4C = ce ACwCr Oran Oo 1 
(b) the liquid concentration on the surface falls to the equil- 
ibrium value immediately at the start of the drying, i.e. 
CameCvdatat -sO candernt=  R< 


i 


With these conditions Newman expressed the solution of equation 


(2.33) as 
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then the expression for the fractional approach to equilibrium, E, 


could be deduced and is, 


co 


Cc -C¢ 
. = Boar t = me oe x iS 4 2 
Co - C; ‘ T2 n=] n2 8 SoS 


2 Z 
gtt/R?) 9 ---------------- (2535) 


Crank [20] has also derived a similar equation for diffusion 
into or from a sphere. 

Calderbank [11] assumed the first term of the series in equa- 
tion (2.35) only to be significant for values of E greater than 0.4 


(long contact times) and obtained, 


E.= 1-6. exp (-n*p_t/R*) | ----~-----—------~~—----=------------ (2.36) 
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where a' is the interfacial area per unit volume of dispersed phase, 


Calderbank derived an approximate equation for web which is 
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and approximated equation (2.37) as, 


-In (1-E) =E=k Gye ----------------------------------+--- (2.40) 


where t is the total residence time of the dispersed phase in the 
phase contacting equipment. 

Geddes [30] also defined k, as in equation (2.37) and assumed 
only the first term of the series expansion in equation (2.35) to be 
Significant. Thus he obtained, 
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where t is the gas contact time. 

Kronig and Brink [51] have considered the effect of internal 
circulation on the extraction of a substance from a spherical drop 
falling in a stagnant fluid, and have derived the following equation 
for the efficiency of mass transfer, by assuming the concentration of 
the diffusing substance, and the resistance to mass transfer in the 


continuous phase to be negligible. 
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where AL Fuh dey Az = 0.73 


r 
r 


1.678, Az = 9.83 
This equation has been accurately represented empirically by 


Calderbank and Korchinski [12] as, 
= aes = 2 2 1/2 
E = [1 - exp (-2.25 Dot t/R°) |" = ---------------------------- (2.43) 


Vermeulen [82] presented an accurate representation of equa- 


tion (2.35) by 


E = 1 - exp (-n’ tea sas SL ee SUNY Son ET PN Re Re So (2.44) 

By comparing equations (2.43) and (2.44) it becomes apparent 
that the Kronig and Brink model of internal circulation results in an 
effective diffusion coefficient of 2.25 times the molecular value. 
Calderbank [11] and Miller [65] report a similar observation. 

Experimental internal drop transfer rates have been found to 
vary from the rigid sphere characteristic for highly viscous liquids 
up to three or four times the Kronig and Brink rate for fluids of low 
viscosity by Calderbank and Korchinski OAE McDowell and Meyers [61], 
Skelland and Wellek [75], Heertjes et. al. [40], and Johnson and 
Hamielec [48 ]. 

Mehta sudesherua [63] studied the gas phase mass transfer 


resistance of a chain of bubbles and found k,a to be proportional to 
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of the liquid seal. 


where w is the gas flow rate and H is the height 
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2.4.2 Studies During Bubble Formation 

Calderbank [9], Thorogood 77ra: and Haselden et. al. [37] 
are believed to be the only ones who have tried to isolate the mass 
transfer into a growing bubble from a continuous phase under gas-phase 
controlled conditions. 

Calderbank [9 ] assumed that as the bubble is formed. from the 
orifice, concentration gradients which may be set up at the interface 
never extend to the centre of the bubble, since this is continuously 
being fed with fresh gas. Thus it was permissible to use the Higbie's 


penetration theory for diffusion into a semi-infinite body, i.e. 


where t' is the mean age of the interface. 

Since there is a high degree of circulation within the bubble 
t' is some small fraction of the total time t from the commencement of 
bubble formation. Calderbank [9] assumed this fraction to be C. t) 
and thus expressed the variation of the gas phase transfer coefficient 


during the life span of the forming bubble as 


where x is a constant characteristic of the system, e.g., x = 1 for 
transfer from a rising gas bubble under continuous phase controlled 
diffusion. 


On this basis Calderbank [9] derived an expression for the 
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maximum efficiency attainable during bubble formation when the bubble 
could be assumed spherical and when bubble volume increased regularly 
with time. The latter condition implied that flow rate of gas to the 


centre of the bubble was constant. The equation is 


oat | CLE vt 


where nes 6 (x Do/m)* net ais the 


[amwy/3 
formation time of bubble, and w is the flow rate of fresh gas to the 
bubble. 
Thorogood [77] has modified Calderbank's above mentioned 
approach for the case when the flow rate of the gas into the bubble 
increases as the bubble grows in size. He assumed the level of turbu- 


lence inside the bubble to remain constant during the bubble growth 


and consequently represented ka during formation as, 


where Dac is effective gas phase diffusivity and te is the time of 
formation. 

He also assumed the height of the top of the bubble above the 
surface of the orifice to be proportional to the time elapsed from the 


commencement of bubble growth and derived the expression for efficiency 


as: 
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2k, + (dh/dt) 
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34 
where h is the height of bubble at any time. 


Thorogood's [77] assumption of the height of the bubble in- 
creasing at a constant rate has been supported experimentally by 
Haselden and Thorogood [37] who studied the distillation of the oxygen- 
nitrogen-argon system on a perforated plate. 

There are indications that a significant amount of mass trans- 
fer may be occurring in the bubble formation regime. Calderbank [9] 
carried out experiments on evaporation rates of water and some organic 
liquids into air wherein air was present as the dispersed phase and 
obtained the mass transfer at zero submergence by extrapolation. The 
amount of transfer was considerable and was interpreted as occurring 
during bubble formation. Haselden and Thorogood [37] calculated the 
fe seat arieaser during bubble formation as the difference in vapour feed 
rate and the computed vapour throughput from the observed bubble sizes 
at detachment and the bubble frequency, and concluded that 20 - 50% of 


the possible transfer may occur during formation of the bubble. 


2.4.3 Mechanisms for Mass Transfer During Drop Formation 


Heertjes and de Nie [39 ] have proposed two models for the 
mechanism of mass transfer to drops during formation when the resist- 
ance in the continuous phase is negligible. Both of these are based 
on no internal circulation in the drop and are applicable when the 
relationship eee area and volume of a drop at any time follows 


the equation: 
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where A(t) and V(t) are the surface area and volume respectively for 
the drop at any time t, A, is surface of the restdrop, and B is a 
proportionality factor independent of the flow rate and time but de- 
pendent on the orifice diameter. 

Also, mass transfer is assumed to take place by unsteady state 
diffusion into a semi-infinite medium. 

The first model assumes mass transfer through a stretching 
surface and results in 
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where E is defined as the efficiency of mass transfer and is equal to 
actual mass transferred divided by the total possible mass transfer. 
The second model assumes mass transfer through a surface ex- 
panding by the formation of fresh elements such that there is no mixing 
between elements of different age and predicts the efficiency of trans- 


fer as: 


A : Dt, 1/2 
E= [2 o + 3B] 1 seatat alr pe A ele eeat ge acire cricket ia tars Sic Wath or wa wm (2.52) 
: 1 


Groothius and Kramers [32] have also derived a result analogous 
to (2.52) but have not considered the effect of the restdrop. 

The experimental results of Heertjes et. al. [39] and Groothius 
et. al. [32], on transfer to isobutanol drops in water and sulphur 
dioxide absorption by water drops respectively, are in fair agreement 


with the values predicted from equation (2.52). 
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CHAPTER 3 


THEORY 


3.1 General Theories of Interphase Mass Transfer 


The various mass transfer mechanisms which have been postulated 
are based on Fick's laws for diffusion processes. Fick's first law is 


applicable for steady state conditions and is, 


Fick's second law for unidirectional diffusion applies to unsteady 


state mass transfer and can be written as, 


Both of these laws are applicable when there is no fluid motion. 


3.1.1 TheeFilimeTheory 

The film theory was first proposed by Lewis and Whitman [56] 
in 1924 and the basic concept involves occurrence of mass transfer into 
or out of a moving fluid phase by molecular diffusion through a thin 
film of stagnant fluid at the phase boundary. Mass transfer through 
the stagnant film is presumed to be considerably slower than in the 
bulk phase and controls the overall transfer rate. Steady state 


mass transfer is assumed to take place and the diffusion is assumed to 
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occur only in a direction perpendicular to the interface. 
With these assumptions Fick's first law holds and equation 


(3.1) on integration yields, 


Defining the mass transfer coefficient as, 


ens) (meme Webcce bos’ —She-rr ese sn (3.4) 


one obtains from a comparison of equations (3.3) and (3.4), 


where, D and 6 are molecular diffusivity and stagnant film thickness 


respectively. 


3.1.2 The Penetration Theory 


The penetration or surface renewal theory was introduced by 
Higbie [41] in 1935. The key difference between the film theory and 
this theory is the applicability of the latter for unsteady state 
transfer processes. The main assumptions are: 
(a) discontinuous contact between phases 
(b) fresh fluid is brought to the phase boundary for each new 
contact and the length of time each new element of fluid 


remains at the interface is constant. 
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(c) the fluid, during the period of contact, remains quiescent, 
and mass transfer occurs by molecular diffusion in a direc- 
tion perpendicular to the interface. 

(d) the fluid field is semi-infinite in depth 

(e) the transfer is entirely by molecular diffusion. 

Fick's second law (equation 3.2) applies and the conditions 
required for solving are: 

(i) the initial concentration of the absorbing material is uni- 


form throughout the fluid; 
G(Z;,0) == C. 


(ii) the concentration of dissolved material deep in the fluid 


remains constant at the initial concentration; 
C(~,t) =C 
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(iii) the concentration of solute is constant at the interface; 


c(O,t) = C, 


The solution of equation (3.2) is, 
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If concentration in equation (3.6) is differentiated with 
respect to Z, and the differential is evaluated at Z = 0, the result, 
when multiplied by molecular diffusivity, D, yields for the instantane- 


ous mass flux for solute crossing the phase boundary, 


D 
We tt) = Oy) (C; Eee maces She geet eee +e Seeeeed—ta_he_oe (3.7) 


When this expression is averaged over time of exposure, which 
Higbie assumed to be a constant for all the fluid elements contacting 


a phase boundary, the result is 


SDR 5s (oul Ce ycbertee teats ee te (3.8) 
The total amount of solute absorbed in time t can be expressed 

as 

fey ne PS CORT 9 2 cect at an hina (3.9) 
Equations (3.8) and (3.9) give, 

kK=2V7(2)  -------------------------------------------------- (3.10) 


where, t is the exposure time of fluid elements. For a bubble, of 


diameter d, rising in a liquid with a velocity U, this exposure time 


is often taken as, 
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3.2 Continuous Phase Controlled Mass Transfer 

To analyse the mass transfer in the continuous phase around a 
growing bubble two models are proposed based on different mechanisms 
for the expansion of interface. In the first model the surface is 
assumed to stretch, that is all elements of surface have the same life 
history, and in the second model the surface is assumed to be composed 
of elements of different ages but no transfer between these elements 


is assumed. 


3.2.1 Mass Transfer Through the Stretching Surface 


This model was initially proposed by Beek and Kramers [3] and 
an attempt is made to clarify their derivation. 

The expanding surface is assumed to grow spherically due to a 
time dependent source at the origin. The origin of the spherical 
cavity remains fixed in space throughout the formation period. The 
only motion of the continuous phase is assumed to be in the radial 
direction which is a consequence of the surface growth. Further, the 
concentration of the solute at the interface is presumed to remain 
constant and far away from the interface the concentration of the same 
component is zero. Any increase in surface is assumed to occur by 
stretching and hence all elements of surface are exposed to the con- 
tinuous phase for the total formation time. The fluid underneath the 
surface is assumed to move with the velocity of the interface. They 
assumed that the normal gradient of the tangential velocity is zero 
at depths where diffusion takes place. 


Let the radius of the growing surface at any time be, 
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Ete), =.(35 a-t 

Pp 
where “e is a point source defined as (wW/47) in which W is the flow 
rate into the bubble. 


Volume of bubble at time t is, 


The strength of the source is given by, 


d 
4 Ta, (p+) ee (3.13) 


The equation for the non-steady diffusion taking place outside 


the surface of the growing cavity is, 


so 22? of) _y Bieat os te "relers ce the disteme of 2 (3.14) 
This equation applies in the region r > R and for p >-l. The 
first term on the right hand side represents the transfer due to molec- 
ular diffusion and the second term represents the convective transfer 
in the radial direction due to the bubble growth. 
we is the radial velocity at r from the centre of bubble and 
can be expressed as, 
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The boundary conditions are, 


p+1. 1/3 


) for C20 


Pape sC=20) abr =°(3 Bas 
i.e. the concentration at the interface remains constant. 

(43). 0. = O-at r= © for t>0 
i.e. the concentration of material being absorbed remains 

zero in the bulk liquid. 

(iii) C= Oat r>O fort = 0 
i.e. the initial concentration of material to be absorbed 

is zero throughout the liquid. 


Introduce the transformation, 


r is time independent as it refers to the distance of a 
point outside the bubble surface from the centre of bubble, which does 
not change with time. R is the bubble radius and is constantly chang- 
ing with time. Hence P becomes time dependent. 

The radial velocity as expressed in equation (3.15) is measured 
from a fixed point, namely the centre of bubble. Once the transform- 
ation (3.16) is introduced the radial velocity becomes zero as it is 


now referred to a time dependent variable. 


| ane 
O78 aaokatbaoe + ana at 
ae 


: 
0 <2 xo ENS ima 4H ED =4 or ’ 2 k 


.?opjeoo> entamsr asornitosnt dd 34 dolsszinssnoes st .a.t 
cn a ; ss 4 
O<3 yo!” = + Js O =D ae. 


ackamet bedteeds goited [sl1stsm to noljszinsomoo 93 ay de 


-htuprl Alud saz at O19s 


O= 3 701 OS 1 In O = JD Cee 


. bedzoads od of Letsastem to nolgextasonon fstrini arly .a.5 


_ i AZ 
-bhupkt ofa tvodguends oveg OE 


ia 


fotinamyotens+ ays ssubhosaad. R 


Aint 
*, ae} 
af ee ee ee ee ee ae ee Le EE EO me ee ee a st tn =: ; 


f 4¢ 
a to sompgeath eda of 2reisy til as insbasqsbal omks eb t 


x ‘ ; re ~~ 
apab dotdw .sididud to sitasd of2 mort eoaltue siddud e823 Shiegae 
a > 


-unstio yitaesano> et bre subbss elddud of} et A .amks 3d 
: ,@ 


.inebrnegsb salt eeaoosd § sonelh . 


bemweaum et (¢€i.f) oolisvpo al boaaaqxo eb xitaotav a 


; em or | pie sed 

—~wyobensi1 ait san0 .alddud to S1Inbo 1). ylomen ,Intog be 
: ' : a 

fees 


va 


at 3i sp otes temooed yIFooloy fetbax sis boouher Sia 


4 


2idslisy staebnsqeb omit 


Ore 

dc dc 3 pt1. 2/3 1 

oT - | e + cee I seston ould Si wor te aes oN oes Wile des ias tach ous ins oae eid obs catia aorta 
ar 0 (op “bus ) 02 


970° 3 pt1.-2/3 ay 3 pt1.-5/3 
aoe 205 (tr hee ) - 2 CE Jact ) 
2(p3 + speed tykhd 2(03 + 3a tPtly4/3 
= a eee eae 
p2 ps 
3 pt1.1/3 3 pt1.4/3 
926 i dc 200 ="+53aF€ ) 1 2(9° + hab: ) 
dr2 90 p2 05 
ke 1 +14 
a age (0° +3 tP ) [3 5 -esererenctsigare anaes ae 
Equation (3.14) now reduces to: 
pt1,1/3 
Pi a hasantes) 
dc sVop oan 3 pt1,1/3 dc 
inioass aE Sa ae + pes 
Tee U2 dp p2 ces ate ed 30” p2 
2(p° + 3a ern j26 3 pt+1.4/3 
fe ae ] + dp2 ° b (Oo i+: 3a t ) 
1ai3 
Ph tga. +) TEAL a 
r 5 Lo 3p2 oa ptl 2) 
p 0? sae t 
or, 
D(p® + 3a pet, 4/3 2 Sane 
eee epee heey: Pe enh eee eo 
a 2 
ot 05 dp Poaiga +P 1 dp 


The new boundary conditions for equation (3.19) are: 
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The solution is greatly simplified by restricting it to a 


region, 


Pye itself is 


which may not be correct at a very early stage as Sri 
very small. But this is relatively unimportant as the surface is 
small and the mass transfer at this stage may not be significant in 


comparison to the rest of the time. 


With this condition equation (3.19) reduces to, 
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Solving this ordinary differential equation with boundary 


conditions (3.20), we get 
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The rate of mass transfer at time t is given by Fick's law 


and is, 
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The total amount of mass transfer during formation can be 


calculated by integrating equation (3.29) over the entire life of the 


bubble 
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The concentration gradient at the interface can be found 


from equations (3.28), (3.23), (3.16) and (3.32) and is 


Equation (3.31) on substitution for R from equation (2.32) 
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For a time independent source at the centre, 
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and from equation (2.32) the radius is, 
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The bubble radius at the end of the formation period is given 


by 
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Substitution from equations (3.36) and (3.38) into (3.35) yields 


for the mass transfer during the entire bubble formation period, 
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Popovich et. al. [67] have summarised the various models 
available for mass transfer during single drop formation and report 


that Ilkovic's model results in the same equation as (3.39). 


The average mass transfer coefficient can be defined as, 
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For spherical growth of the bubble, 
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Substitution for = from (3.38) into above gives 
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Equations (3.39), (3.40) and (3.44) can be simplified to give, 


This close resemblance with the general solution of the 
penetration theory is not surprising as the mean exposure time of the 
elements was assumed to be te in this derivation. 

Equation (3.39) is valid only when the depth of penetration 


is small compared to the radius of the bubble as it was one of the 


assumptions (3.21) in the derivation. Quantitatively when, 
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where, (0/2) is assumed to be the effective film thickness. 
Equations (3.47) and (3.48) can be combined to give the final 


condition for the application of equation (3.39) as, 


p3/2 1/2 o 


3 
Re 
The mass transfer coefficients during bubble formation and 


bubble rise can be compared from equations (2.22), (3.11), (3.45), and 


(2.29). The result obtained is, 


ka (formation) d 
ke (rise) meh GG 


which implies that the mass transfer coefficient during bubble form- 
ation is always lower than its value during bubble rise as (d,/2%) is 
always less than one if there is no bubble coalescence at the orifice. 


The limitations of this model are as follows: 
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(a) the effect of continuous phase velocity is not considered 
(b) spherical bubble growth is assumed which is the case of 
minimum transfer as any deviation from sphericity would in- 


crease the area and hence the transfer. 


3.2.2 Mass Transfer Through An Interface Expanding by the Formation 


of Fresh Elements 

In the previous model each element of the interface was assumed 
to have a life span equivalent to the total formation time. 

Calderbank and Patra [15] have represented the growth phenom- 
enon by assuming that the additional interface formed in the course 
of the time is entirely fresh and that there is no exchange of matter 
between elements of different age. This shall be used as the basis for 
our model. All the other basic conditions of Higbie's penetration 
theory are assumed to apply. 

The general solution of the penetration theory gives the 


instantaneous rate of mass transfer per unit area as, 
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By integrating, the amount transferred in time t per unit area 
is 
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The total amount of mass transfer from the growing bubble 
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51 
starting from an initial area S(O) consists of: 
(a) transfer from the surface elements having a life at the inter- 
face equal to the total bubble formation time, te. 
(b) transfer from the surface elements having a life span of 
(t, - 0); where 9>0. 
Since no surface element is assumed to leave the interface once 
it reaches there, it is imperative that area S(0) will have a total age 


of ty and the amount transferred from these elements is obtained from 


equation (3.51) as, 
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Any subsequent surface element added to the interface during 
bubble growth has a unique life span uncommon to other elements. If 
one such element is added at time 0>0, the amount transferred from it 


per unit area from equation (3.51) will be, 
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The contribution to the total mass transfer from all these 


unique elements can be expressed as, 


The total mass transfer during bubble growth is found by adding 


equations (3.52) and (3.54) and is, 
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Substituting from equation (3.53), we get 
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Assume the area of the growing bubble to vary as, 
rn na eee (sea) 
where, S(9) = Surface area of bubble at any time, cm’. 

Sy = Surface area of rest bubble remaining at the orifice 


after bubble detachment, cm”, 
Re = Slope of area - time curve, cm/sec. 
With this substitution equation (3.56) gives, 
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where, Dec) = Total mass transfer during bubble formation, gm-moles 
Cc, = Concentration of material being absorbed at the inter- 
face, gm-moles/cm’. 
G. = Concentration of material being absorbed in the bulk 
liquid, gm-moles/cm?, 
dD, = Liquid phase diffusivity, cm7/sec. 
te = Time of bubble formation, sec. 
Sy = Surface area of rest bubble, cm*, 
K' = Slope of area-time curve, cm’, 


The main limitations of this model are that the effect of 
liquid velocity on mass transfer is not considered and the penetration 
theory is assumed to apply. 


Also, if ke is defined by 
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then substitution from equations (3.57) and (3.58) gives the average 


mass transfer coefficient during bubble formation as, 
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where, ke = average continuous phase mass transfer coefficient during 
bubble formation, cm./sec. 


Calderbank and Patra [15] had also assumed a similar model and 


had derived an expression for Ke by defining it as, 


mm, = ke AL te (Cc, - C_) 0 -eeenneee- eee = = ++ (25.26) 
which gave, 
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They had proposed an empirical correlation for the variation 


of bubble surface area with time as 
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which implies that the area at time zero is equal to zero. If one 


neglects S_ in equation (3.60), the solution obtained is, 
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Equations (2.27) and (3.61) are identical except for the value 
of the constant. In Calderbank's equation (2.27), the value of kK is 
low because of the way he defined ke in equation (2.26). 

To have a true comparison it is imperative to define ke in the 


same fashion. Thus, by changing equation (2.26) to 
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equations (3.59) and (3.62) become compatible. 
Substituting for A in equation (3.62) from equation (2.9), 


integrating and comparing the resulting equation for m, with equation 


f 
(2.25) one finds that Calderbank's solution is expressed by the fol- 


lowing equation instead of equation (2.27) 


It can readily be seen that equation (3.63) is very close to 
the solution derived here, equation (3.61). This difference in the 
values of constant can be due to different area-time relationships 
used in the two solutions, equations (2.9) and (3.57). 

From equations (2.22), (3.11), (3.60), and (2.29) it can be 


shown that, 
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which for very low frequencies of bubble formation predicts the mass 
transfer coefficient during formation to be smaller than its value 


during bubble rise. 
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CHAPTER 4 


EXPERIMENTAL 
4.1 Choice of Systems 

The carbon dioxide-distilled and dimineralized water system 
was selected to determine the liquid phase mass transfer coefficient 
during bubble formation. The choice was justified as, 

(i) the solubility of carbon dioxide in water under atmospheric 
conditions is small and so the resistance to mass transfer 
can be assumed to be in the liquid phase; 

(ii) the diffusivity of carbon dioxide in water can be assumed 
constant at a particular temperature; 

(iii) the heat effects are negligible as the heat of solution of 
carbon dioxide is small. 

The ammonia-nitrogen~-distilled dimineralized water system 
was chosen for a study of mass transfer in the dispersed phase for 
the following two reasons: 

(i) the major mass transfer resistance can be assumed to be in the 
dispersed phase as the solubility of ammonia in water is very 
high; 

(ii) the nitrogen provides transfer resistance in the dispersed 


phase and its transfer to air saturated water can be neglected. 


4.2 Basis of the Experimental Method 


In general, the mass transferred during the formation of a 
bubble was determined as a difference between the gas flow rate to the 


orifice and the product of the frequency of bubble formation and the 
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mean volume of the detached bubbles. Mathematically 


Mass transfer =W - n Ve wr (4.1) 
where, W is the gas flow rate in cc/sec 
n is the bubble frequency in sec 


Ve is the mean volume of the detached bubbles in c.c. 


The validity of the experimental method employed was checked by 
using the air-water system wherein no mass transfer should occur if 


water saturated air and air saturated water are used. 


4.3 Experimental Set-up 


4.3.1 Carbon Dioxide-Water System 

A schematic flowsheet of the experimental set up is shown in 
Figure 2. Bone dry carbon dioxide gas (minimum purity 99.8%) from a 
high pressure cylinder (1) was saturated with distilled water in a 3 
inch diameter lucite absorber (2) having a 3 feet packed height of glass 
beads (1/4 inch o.d.) and was then metered through a soap film meter 
(3). The gas pressure at the soap film meter was regulated by a press- 
ure regulating valve on the gas cylinder and fine control was subsequ- 
ently achieved by a diaphragm valve just before the soap film meter. 
The pressure reading was taken on a mercury manometer (5). The gas 
flow rate was accurately controlled by a micrometer needle valve (4). 


The gas then passed through an adjustable chamber volume (6) to the 
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orifice (7), made of 316 stainless steel tubing 0.408 cm. inside 
diameter. A teflon cover was provided at the top of the orifice to 
make it non-wettable with water. The orifice was installed in a lucite 
bubble column (8) having an internal diameter of 20.2 cm., a wall 
thickness of 0.635 cm., and a total height of 100 cms. 

The water was stored in a lucite reservoir tank (11) of about 
100 gallons capacity and was pumped by a centrifugal pump (10) through 
a set of rotameters (9) into the column (8). The flow of water and gas 
through the column was co-current. The water overflow from the column 
(8) went back into the reservoir tank (11). Polyvinyl Chloride was 
used as the material of construction for all the liquid pipe lines. 

Temperature control to + 0.05°C was achieved with the use of 
a Beckman thermometer (13) installed in the water reservoir (11). A 
control circuit was automatically activated by the thermometer when the 
water temperature exceeded the desired value. When this happened a 
small capacity centrifugal pump (14) started and circulated cooling 
water from a refrigerated storage (15)(the temperature in which was 
maintained at about 4°C) to a cooling coil (12) immersed in the water 
reservoir (11). 

The experimental set-up described above is essentially the 
same as that used by Kalra and described in detail in his thesis _49]. 
Some modifications were, however, made towards perfecting the exper- 
imental technique. The most effective change was the introduction of 
a diaphragm valve, for controlling the gas flow rate, just before the 
soap film meter. A very steady gas flow could be obtained even at ex- 


tremely low frequencies of bubble formation (one bubble per two sec- 
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onds). Also, since the gas pressure reading is very important in 


subsequent calculations it was read on a mercury manometer. 


4.3.2 Ammonia-Nitrogen-Water System 

The flow sheet of the equipment used for studying the 
ammonia-nitrogen-water system is shown schematically in Figure 3. 

Ammonia gas from an anhydrous ammonia cylinder (1) was metered 
through a previously calibrated Matheson gas flow rotameter (2). A 
soap film meter was tried to measure the ammonia gas flow but was ab- 
andoned as the soap film kept breaking due to the high solubility of 
ammonia in soap-water solution. Nitrogen gas from a high pressure nit- 
rogen cylinder (16) was mixed with ammonia gas in a T-mixer after met- 
ering through a soap film meter (17). A fine flow control for both 
the ammonia and the nitrogen was achieved individually by installing 
a diaphragm valve in each gas line before the gases enter their re- 
spective flowmeters. The flow of ammonia-nitrogen gas mixture was 
then measured in another soap film meter (3) to allow a check on the 
material balance. The remaining flow scheme is the same as described 


above for the carbon dioxide-water system. 


4.4 The Photographic Technique 

A high speed camera was used to take motion pictures of the 
bubble during its growth from the orifice. A key factor in obtaining 
Heit defined periphery of the bubble on the pictures was the technique 


used for lighting. The photographic set-up used is shown schematically 


in figure 4 in which sketches a, b and c show the layout of the side 
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(A) SIDE-LIGHTS (B) BACK-LIGHTS 
AND CAMERA 


(c) PLAN OF PHOTOGRAPHIC 
SET-UP. 


LEGEND 

1,2,3. PAPER LIGHT SHIELDS 
4 GROWING BUBBLE 

5,6,8 LAMPS 

7 ORIFICE 

9 CORRECTION LENS 

10 HIGH SPEED CAMERA 

i" CAMERA STAND 

12 CAMERA TABLE 

13 LUCITE BUBBLE COLUMN 


FIG. 4. PHOTOGRAPHIC ARRANGEMENT. 
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lights, the arrangement of the back-light and the camera, and a plan 
view of the entire photographic set-up, respectively. 

The lighting was provided by two side lights (5, 6) and a back- 
light (8) each with an illuminating intensity of 300 watts. The light 
from the former came at an angle to the orifice, whereas it came horiz- 
ontally from the latter. Surplus light was reflected back by aluminum 
foil sheets (1, 2, 3). The inside surface of these foils was black 
which helped in cutting out the light reflections from the lucite sur- 
face of the column to the bubble. 

A correction lens (9) was used in front of the bubble column 
to avoid distortion of the bubble due to curvature of the column. 

This lens consisted essentially of a hollow lucite casing filled with 
water and had the same curvature as that of the bubble column on one 
side and a plane edge towards the camera. A detailed description of 
this correction lens together with the ray diagram is given in Kalra's 
thesis [49 ]. 

A Hycam, model K2004E, 16 mm. high speed motion picture camera 
was used for taking the pictures and it had a maximum film capacity of 
400 feet and a speed range of zero to 11,000 frames per second. The 
camera (10) was supported on a camera stand (11) which had provisions, 
for horizontal movement on a camera table (12), for horizontal rotation 
of the camera, and for adjustment of the inclination of the camera lens 
axis to the horizontal plane. The camera table (12) was fitted with 
wheels to enable the movement of the camera unit as a whole. 

The Eastman Tri-X Reversal film type 7278 was used in the 


camera and was developed as a reversal film meaning thereby that the 
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64 
dark and bright spots of the object appeared as dark and bright spots 
respectively on the film. 

A Milli-Mite timing light generator model TLG-3 was connected 
to the camera and could be used to determine the exact speed of the 
film as the light generator left a mark on the edge of the film at 
adjustable frequencies of 10, 100 and 1,000 per second. 

During the development stages of the above mentioned photo- 
graphic set-up some problems were encountered. Since a circular column 
was used, light reflections on the periphery of the bubble image on the 
film were noticed. However, these reflections were eliminated by using 
the balanced lighting approach illustrated in figure 4. Also, some 
difficulty was encountered in obtaining a good colour contrast between 
the object and background as both gas bubble and water were colourless. 
The solution was found after trial and error and it involved finding 
the right combination of the intensity of lighting, the speed of the 
camera in terms of frames of film advanced per second, the focal length 
of the camera lens, and the aperture of the lens. Very good pictures 
were obtained when, the orifice was illuminated by three 300 watt lamps, 
the camera speed was 600 frames per second, the camera lens had a 


100 mm. extension, and the lens opening was f/11.0. 


4.5 The Experimental Procedure 

Since the experiments involve the measurement of very small 
quantities, a very careful experimental technique is neeced to obtain 
reliable results. The experimental procedure which was followed 


closely during the course of this investigation is described below. 
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65 
4.5.1 The Preparation for Runs 
The equipment was carfully cleaned prior to making any exper- 
imental measurements and several checks were made before each exper- 


imental run to ensure a proper working of the set-up. 


(a) The Cleaning of Equipment 

The cleaning of equipment is essential as any contamination in 
the system could conceivably affect the reliability of results. 

After the initial assembly of the equipment, the system was 
left running with air bubbles being formed in circulating water for a 
few hours. The water was then drained completely and the procedure 
repeated with a fresh batch of distilled water. After repeating this 
for three days it was assumed that all the suspended impurities like 
PVC chips etc., were completely removed. At this point distilled 
water, containing detergent, was introduced in the water reservoir and 
recirculated through the water pump, rotameters and column for a day. 
This was done to remove any traces of oil entrapped in the system. The 
equipment was rinsed free of detergent by recirculating fresh batches 
of distilled water about fifteen times. 

The column was then opened and glass beads, previously cleaned 
by immersing in concentrated chromic acid for a couple of days and 
rinsed thoroughly with distilled water and dried, were introduced into 
the bottom section of the column. The purpose of these glass beads is 
to obtain a flat velocity profile for the liquid in the column. The 
teflon orifice cover was ain cleaned inside and outside by detergent, 


washed with distilled water and dried. It was then attached to the 
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orifice tip and the column was completely sealed. One fresh batch of 
dimineralized distilled water was circulated through the system and 
drained out. 

The equipment was then considered ready for a calibration run. 

After every run the water was drained to below the level of 
teflon orifice cover in the column. Before any subsequent run the 
water reservoir was completely drained and the system flushed with a 
fresh batch of dimineralized distilled water. The flush water was dis- 


carded completely and the system was refilled with dimineralized water. 


(b) Checks Before Any Run 
Before proceeding with any experimental run the following 
checks were made: 

1. The system was checked for leaks in the gas lines by closing 
the needle valve (4) in Figures 2 and 3 and bringing the gas 
lines to a pressure of about 5 PSIG by introducing air and then 
leaving the isolated system under pressure overnight. 

2. The quality of dimineralized distilled water being fed to the 
water reservoir was continually checked. The satisfactory 
limit was set at a maximum salt concentration of 0.05 ppm. The 
cartridges in the Bantam Dimineralizer were replaced if this 
limit was exceeded. 

3. The camera was checked to make sure that it was, in working 
order, and running at the indicated speed on the dial. 


4. The gas flow rate was checked by taking measurements on the 
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soap film meter to ensure that a steady flow could be obtained 


at the desired frequency of bubble formation. 


4.5.2 The Calibration Run With Air-Water 

The supply air pressure was regulated to the desired value and 
the diaphragm valve (16), Figure 2, was set to give a steady pressure 
of about 2 PSIG on the discharge end. The air was fed to the saturator 
(2) slowly to prevent a surge of water through the gas lines. The mic- 
rometer needle valve (4) was left wide open during the initial stages. 

The column (8) was then filled with water at a very slow rate 
and then the rate was adjusted to the desired value. The set point 
of the Beckman thermometer (13) was fixed at 25°C and the temperature 
control circuit was turned on. The flow of cooling water through the 
cooling coil (12) was adjusted manually till the temperature could be 
controlled to within a limit of + 0.05°C. 

The equipment for measuring the bubble frequency - namely the 
laser, photocell, universal counter and the oscilloscope - were turned 
on. The description of the circuit in detail is given in Kalra's 
thesis [49]. The principle of operation is that once the laser beam 
is adjusted such that it intercepts the bubble after detachment, the 
photocell interprets every interruption of the beam as a discrete 
signal and transmits it to the universal counter. The universal 
counter displays the time elapsed between each interruption and this 
can be taken as the time period of bubble formation. This time period 
was brought close to the desired value by adjusting the needle valve 


(4) opening. 
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Enough time was allowed for to attain steady state. The steady 
State in this context is defined as the state of the system when the 
bubble frequency, gas flow rate through the soap film meter (3), water 
reservoir temperature and water flow rate have attained constant values. 

All the variables were recorded and bubbles photographed by 
using the procedure outlined in (4.5.3). The variables were once again 
recorded to get an average value of the variables during the time 
pictures were being taken. 

The water pump, cooling circuit, frequency measuring circuit 
were then turned off. Water was drained from the column to the bottom 
of the teflon cover on the orifice. The gas supply valve was then shut 


Off 


The photographic arrangement was previously described in 
article 4.4 and shown schematically in Figure 4. The 'loaded' camera 
was adjusted to such a position that the lens axis was exactly horiz- 
ontal, and perpendicular to the correction lens. This is necessary to 
avoid any distortion of the bubble image. Focussing was done on the 
orifice with back light turned on and the lens aperture fully open. 
The appropriate values were set on the speed control of the camera and 
frequency output of the Milli-Mite cycle timer. The camera was locked 
in this set position to avoid any movement during the run and pictures 
were taken with all the lights turned on after all the variables pert- 


inent to the experiment had achieved stable values. 
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4.5.4 Mass Transfer Runs With Carbon Dioxide-Water System 

The procedure was essentially the same as already described for 
the air-water system in article 4.5.2. The exception was that much 
more time was given to flush the gas lines with carbon dioxide gas 
before filling the column with the liquid. The purpose of doing this 
was to remove any entrapped air from the lines. 

The continuous phase was analysed for carbon dioxide to deter- 
mine the driving force for mass transfer. The analysis of CO2 in water 
was done by back titrations. A blank titration of Be(OH\S 261 GEL on 
against standard HCl solution was done for each experimental run in an 
atmosphere of nitrogen. To titrate the water sample, a known volume of 
it was added to a known volume of Ba(OH)2 solution and the mixture was 
titrated with HCl to a green end point of Bromo Thymol Blue indicator. 
A nitrogen atmosphere was used throughout the titration. The concent- 
ration of CO. in water could then be determined by calculation as in- 
dicated in Appendix C of Kalra's thesis [49]. The water sample was 


analysed both before and after taking the pictures. 


4.5.5 Mass Transfer Run with Ammonia-Nitrogen-Water System 

A constant composition gas mixture was obtained by accurately 
controlling the flow rates of both nitrogen and ammonia gas separately 
and then mixing them in a T-mixer (see Figure 3). The pressures of 
nitrogen and ammonia lines were adjusted to the same value by setting 
the openings of diaphragm valves (18) and (19) respectively. The rest 
of the experimental procedure followed was the same as that for the air- 
water system described in 4.5.2. The time needed to achieve the steady- 


state was much longer in this case as the composition of the gas mixture 
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70 
varied slightly with any change in the setting of the needle valve (4), 
needed for adjusting the frequency of bubble formation. 

The bulk phase, water, was analysed for dissolved ammonia both 
before and after taking the pictures by titrating the water sample 
against a standard solution of HCl. 

Nitrogen and ammonia gas flow rates were recorded individually 
to check the material balance against the recorded gas flow rate of the 


mixture in the soap film meter (3). 


4.5.6 Sample Data Sheets 

Since it is necessary to record all the variables precisely, 
the format in which the data was taken is included for reference. 
Tables 2, 3 and 4 are for air-water, carbon dioxide-water, and 
ammonia-nitrogen-water respectively. These tables refer to three 
actual runs and will be used later to illustrate the procedure for 
calculations. A number of readings were taken for the time period of 
bubble formation both before and after taking the pictures as these 
readings were used to statistically determine the number of bubbles 


which had to be analysed for volume calculations. 


4.6 Analysis of Results 

Step l The measured gas flow rate at the soap film meter was 
Ronvarekd to the orifice conditions of temperature and pressure 
by assuming the ideal gas law to apply. 

Step 2 The arithmetic mean time period of bubble formation 


was calculated from the recorded readings of time period. 
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The number of detached bubbles that had to be analysed 
for volume was determined statistically from the recorded 


readings of time period by using the formula, 
n 
Zo (x, - x)? --ae n-ne + +--+ +--+ (4.2) 


where, N is the number of bubbles to be analysed 
n is the number of data points 
x is the arithmetic mean time period 

x, is the ree time period 

This formula gives a confidence limit of 99%. The 
calculations were aided by a computer program written for the 
purpose by Kalra and documented in his thesis [49]. 

A minimum of ten bubbles were analysed. 

The films were screened to see that they complied with 
the desired specifications of, sharpness of the bubble and 
orifice edges, and contrast between the bubble and its back- 
ground. 

A number of detached bubbles, as determined in step 3, 
were analysed by using a digitizer set-up available in the 
Department of Chemical Engineering, University of Alberta. 
The film was projected on a horizontal screen calibrated in 
two aireccions (X and Y) in terms of standard voltages. It 
meant that any point on the screen could be represented 


discretely by two voltages corresponding to its X and Y coor- 
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dinates respectively. These voltages could be punched auto- 


matically on to computer cards by the use of the digitizer. The 
orifice image on the film served as a reference object thus 
enabling one to convert any point, read on the projected film 
and punched on the computer card, to its actual coordinates 
as follows: 
Say, 10 cm. on screen = 2 volts 

Voltage recorded by digitizer across 

the projected orifice image = V, volts 

Known orifice diameter = D, cm. 
.. Conversion scale for all digitizer readings = a em/volt 
Let, voltage corresponding to X and Y coordinates of a point 


on the film be V, and V2 volts respectively. 


uf 
Then, the X — coordinate of that point = a x Vi> cm. 
D 
and, the Y - coordinate of that point = vy * V2, cm. 


A chosen number of points on the projected image of the 
bubble periphery were read and transferred on to the computer 
cards. A program previously written by Kalra [49] was then 
used to obtain the surface area and the volume of the bubble. 
His program numerically integrates the following two equations 


for surface area and volume respectively. 


2n f £(v) V[1 + (£'(v))?] dv ----------------------------- (4.3) 


wS[£ (v)]*? dv ------------------------------------------- (4.4) 
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Where f(v) is a numerically determined polynomial function 
describing the periphery of the bubble and f'(v) is the derivat- 
ive of f(v) with respect to v. 
The errors for surface area and volume, calculated from 
the program, were stated by Kalra to be within 0.1% and 1.0% 
respectively. 

Step 6 The amount of mass transferred during bubble formation was 
calculated from equation 4.1. 

Step 7 For the mass transfer runs the bubble having a detached 
volume closest to the mean value of all analysed detached 
bubbles from Step 5, was analysed at various points during its 
formation by following the same procedure as in Step 5. The 
areas of the growing bubble as obtained from the computer pro- 
gram were plotted as a function of time and the area under this 
curve was graphically integrated to give a value for vf A dt 

Step 8 The average mass transfer coefficient during bubble form- 


ation was calculated from the following equations, 


a Sa a a (4.5) 
i 
(fomenraty ic, 
t 
Kec * es oe aM i OS pre oee a ae (4.6) 
rf 
en 


Where ke and K are the liquid and overall gas phase 


GC 
transfer coefficients respectively.  ¢ is average mass transfer 


rate as calculated from equation 4.1, and Cc. and Yn are the 


concentration driving forces for the liquid and gas phase con- 
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trolled systems respectively since from the experimental ana- 
lysis during the mass transfer runs it was observed that the 
concentration of solute in the bulk phase was negligible. 

For the COz-water system Cc, was taken as the solubility 
of CO2 in water, and for the NH3-N2-water system y, was 


defined as, 


Where Ys is the mole fraction of ammonia in the gas feed to the 
orifice and Ye is the ammmonia concentration in the detached 
bubble and is calculated from material balance as shown in 


appendix A.4. 


4.7 Precautions 


A number of precautions were observed during the course of this 


investigation and are strongly recommended to obtain reliable results. 


(a) 


(b) 


The chamber volume should be held constant for all the exper- 
imental runs and the criterian to select this constant value 
should be to obtain a stable bubble frequency. 

The centrifugal pump for water circulation should always be run 
at a certain discharge head. Failure to do this could: result 
in air bubbles being entrained and circulated through the 
column. To ensure that the pump runs properly the valve on the 


by-pass line to the water reservoir from the pump discharge was 
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kept partially closed. 

While draining water out of the bubble column, care should be 
taken to make certain that the level in the column does not 
fall below the top of the glass beads. This prevents air from 
getting entrapped between the beads. 

The purity of the water used should be consistent for all the 
experimental runs. 

During the experimental run the water temperature should be 
controlled very accurately. 

All titrations should be done under a nitrogen atmosphere. 
About half an hour's warming up time should be allowed for 

the universal counter. This corresponds to the minimum time 
recommended by the manufacturer. 

Since the barometer used for recording the atmospheric pressure 
was in a different laboratory, the exhaust fans in both the 
laboratories were switched off during the experiment. 

The level of soap solution in the soap film meters should be 
as high as practicable. This avoids undue pressure fluct- 
uations in the gas lines when the soap film is being formed. 
Before proceeding with the run a quick visual observation 
should be made of the system to ensure that there are no en- 
trapped air bubbles in liquid lines or rotameters and that 

the bubble growth is steady. 

Enough time should be allowed for the system to stabilise after 


any minor change is made. The system response is very slow 


nls e ;" <M P= x 7 
»tddad silg te Jub. aww Brdseher sty 


d bios Jay , Glo: 
| “aon 
} th. amdifoy sional tavel af? 2609 «telaeo wales a0 now " 
= 
id > la 
i iT iq alii Shuasd. & ‘afta dy qos ods wot od it La 
Pc > 
-o4 af) osewiad begentine gal loeg 
Ys » 7 van) 
“4. le 40% anetabendp od biuode boa 1adBy wns jo, yuri wen 
N . a Nd 
aciturt Ladies 7 
: i ) w9aIkW ma: ‘La jesniyeges yy Be pnt 
i ‘ i j 3 
-VileJesIUOPE Tey vel Lenses 
; n7a fn ta! : bo anol ad bipode enotder3id ih 
alle sd hlanide t 1 Zoos 'sy0H op thed sain 
a OF ; a . i & . 
' r Ta - eT Ita ou a - , ALND ime vw au alia 
a oe ~ 
. = : , yi 
12739247 RaN 809 Fo PRDHaReEE 
* Ser “~ ty rv “ ‘ ry tras ‘ wei | eel aoolk : rf 
> * Pa §% . ' a ‘ 
- ; } os 
sii} 30 n 6&1 IG yA {f unsts rtib & al ~ iY 
4 rf 
neath qxo att wutoe | hotptiwe s¥ew Srl sede 
cae b ( q 7: Teh | ni cs) 4; aac a 2a Leann ad) 
g ; oy: 
¢ : gf j : o , ~_ 
~35ut% : ih @biteve rit -o1 33 jJ0a7%g =@s tio bal. a 
DIT yl6sy aoa silt uety weolhl tsp Sis at Bomhaee 
4 ey yy 


* > 
netjeveeedo Lr 


We by) Buse 97S ind eodtens a3 wiouye aid- ip ene 


$405 bas 


wore heade: sinad a muiert ‘was aa bawe iis 


¥) 


a 


ube 


: 


rely Wo bi & aT eval as ew - eeiaieidad Bad 


«ta tener 30% 


Wee: eee al i jenngabt 


idee 


7.7 


an 


te sont tuptt al eased 3 


wl his da a asain 


i P 


2) oe. 
id raed om, 
b le 


pete nt 
y 4 a 


nee 


ti i- 


(1) 


(m) 


(n) 


(o) 


(p) 


(q) 


(r) 


80 
because of the low gas flow rates under study. 
It is absolutely necessary to align the camera properly as 
described in the photographic procedure. Any slight inclinat- 
ion of the lens axis to horizontal could cause a distorted image 
of the bubble. More so, because no reference source was used 
in the vertical direction. 
The camera and table should be as free of vibrations as 
possible. 
There should be no gap between the correction lens and the 
outer surface of the bubble column and the water in the correct- 
ion lens should be free of entrapped air bubbles. 
The camera lens should be focussed on the orifice edge. 
The camera lights should be turned on for as little a time as 
possible in order to reduce heat effects. 
While analysing the films, the digitizer calibration should be 
checked frequently in order to detect any drift which might 
have occured. 
When projecting the film on the screen, the projector should 
be exactly horizontal to the ground and the projector lens 
axis should be at 45° to the mirror surface. This is necessary 


to avoid distortion of the bubble image. 
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CHAPTER 5 


EXPERIMENTAL RESULTS 


oye l General 


(i) 
(ii) 


(iii) 


runs: 
(i) 
(ii) 


(iii) 


The results of this investigation are subdivided into: 
Calibration runs with air-water 

Mass transfer runs with carbon dioxide-water 

Mass transfer runs with ammonia-nitrogen-water 


The following parameters were the same for all the experimental 


Orifice inside diameter = 0.408 cm. 
Orifice outside diameter = 0.738 cm. 
Water height in the bubble column above the orifice 


tip = 91-..7 cm. 


5.2 Calibration Runs with Air-Water 


The runs with air-water were necessary to check the validity 


and accuracy of the experimental technique used and a number of them 


were made over the entire duration of the test program. The first 


calibration run was done before starting the mass transfer runs with 


carbon dioxide-water, the second run was done about half way through 


COz2-water runs, the third run was completed after the technique 


for ammonia-nitrogen-water runs was perfected, and the last run was 


done after finishing the experiments. 


The results together with the experimental conditions are 


presented in Table 5. Table 2 is a detailed data sheet for one of 
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these runs (A2). 


These results indicate that the bubble volume can be deter- 


mined with an accuracy of + 2%. 


5.3 Mass Transfer Runs with Carbon Dioxide-Water 

Initially experiments were conducted with carbon dioxide as the 
gas phase and water as the continuous phase for gas bubble frequencies 
of 3 and 5 bubbles per second. Four different water flow rates vary- 
ing from 0.34 cm./sec. to 1.09 cm./sec. were studied. These runs in- 
dicated the percentage mass transfer occurring during bubble formation 
to be in the range of 0.5% to 2.5% which fell within the experimental 
uncertainty of + 2%. The results are tabulated in Table 6 and are for 
qualitative use only. 

On this basis it was decided to study the mass transfer for 
bubble frequencies of 1 and 0.5 bubble/sec. and water velocities of 
0.3 to 2.5 cm./sec. The choice of lower limit of bubble frequency was 
dictated by the film capacity of the high speed camera. The camera 
could not take more than 400 ft. of film which could only record about 
12 bubbles at 1/2 bubble per second and 600 frames/sec. camera speed. 
By statistical analysis it was repeatedly found that approximately 10 
bubbles had to be analyzed. The camera speed could not be reduced 
below about 600 frames/sec. as it could not record the bubble at de- 
tachment at speed below this value. 

The upper limit of water velocity was set at 2.5 cm./sec. as 
the bubble column started overflowing above this velocity. The lower 


limit of 0.3 cm./sec. was chosen as it was thought to be the minimum 
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water flow needed to maintain a constant water temperature in the 
column and reservoir. 

Table 7 gives the main variables studied in each run. The 
data are recorded in Table 8 and the final results are presented in 
Table 9. Table 3 is a detailed data sheet for run number M8 and sample 
calculations for this run are presented in the appendix (A.3). 

The surface area and volume of the detached bubbles as obtained 
from the films are listed in Table 10 and Table 11 respectively. They 
give an indication of the digitizer consistency. 

The surface area of the growing bubble is presented in the 
appendix A.5 (Tables 18 to 24) as a function of time and was found to 
increase linearly with time during most of the formation period. The 
main deviation from linearity was observed towards the end of the form- 
ation period. This region corresponds to the 'necking' period. The 
bubble becomes elongated and forms a neck just before detaching from 
the orifice. A portion of the forming bubble remains at the orifice 
and the rest moves upward through the column. The former is called 
the 'rest bubble' and the latter the detached bubble. Two typical 
area-time curves are plotted in Figure 5 which correspond to runs M8 
and Mll for 1 bubble/sec. and half bubble/sec. respectively. The 
area of the rest bubble was added to the area of the detached bubble 
since at the instant of bubble detachment both the rest bubble and the 
detached bubble are contributing to the mass transfer. Quantitatively 


this relationship can be approximated by a linear function of the form, 
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*These runs were used for deciding the operating range only 


Table 7: Range of Variables Studied for CO2-Water 
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—:— EQUATION 5.4 FOR RUN M 8 
o——e EXPERIMENTAL POINTS FOR RUN MII 
——-—- EQUATION 5.4 FOR RUN MII 
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2 


where, A = Surface area of the growing bubble at time 't', cm 
A, = Surface area of the bubble at time zero, cm? 
k' = Slope of area-time curve, cm*/sec. 
t = Time at any instant of bubble formation, sec. 


The criterion for selecting the value of slope of area-time 
curve was that fAdt should be the same as that determined exper- 
imentally. 

The slope k' varied from 0.62 to 0.75 cm*/sec. with a mean 
value of about 0.66 cm*/sec. for runs with a bubble frequency of one 
per second. For a frequency of half bubble per second the corres- 
ponding values of k' were from 0.31 to 0.37 cm*/sec. with a mean value 
of about 0.33 cm*/sec. A new constant in equation (5.1) could be used 


to account for the effect of gas flow rate and this gave, 


ee ae NU eee oie ie ae cee > (5.2) 


where, 'n' is the frequency of bubble formation oneal 


AS in equations (5.1) and (5.2) could be approximated as, 


Pee ee. See WO LUS We 22 Se Sees 2 tee see Ceee ese (389) 
fe) 4 “o 
where D is the inside orifice diameter. Hence equation (5.2) 
fo) 
becomes , 
A=tp2+0.66 nt  ------------------------------------------ (5.4) 
Nagi 03 


This equation is plotted in Figure 5 for comparison with the 
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experimental values. As mentioned before equation (5.4) is not valid 
for the period just prior to bubble break-up. 

The volume of the growing bubble was also found to vary line- 
arly with time. The data are presented in appendix A.5 (Tables 18 
through 24) and typical curves are illustrated in Figure 6. The only 
exception to the linear relationship was in the initial 10% of the life 
of the bubble. In this period, from visual observation of the film, 
the bubble was oscillating considerably. The volume of the rest bubble 
was again added to the volume of the detached bubble. 

The quantitative relationship for the volume of the growing 


bubble can be approximately written as, 


where, V = Volume of growing bubble at time 't', cm 


ke Slope of volume-time curve, cm?/sec. 
t = Time at any instant during bubble formation, sec. 
The slope k' varied from 0.085 to 0.096 cm*/sec. with a mean 
value of about 0.091 cm*/sec. for a bubble frequency of around one 
bubble/sec. And the value was 0.0425 cm?/sec. for the bubble frequency 
of about half bubble/sec. 


Equation (5.3) could be generalized to account for the effect 


of frequency and the resulting form is, 


-1 
where, 'n' is the frequency of bubble formation, (sec ) 
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The liquid phase mass transfer coefficient was found to be 
nearly independent of the gas flow rate and the liquid velocity. How- 
ever, as expected, the percentage mass transfer was observed to be 
inversely proportional to the gas flow rate. These results are tab- 


ulated in Table 9. 


2:4 Mass Transfer Runs with Ammonia-Nitrogen-Water 

The rate of mass transfer from ammonia-nitrogen bubbles forming 
in a continuous water phase was measured experimentally. The effect of 
gas flow rate, or the bubble frequency, on the coefficient was studied 
while the water velocity in the column was maintained constant at 
1 cm./sec. The gas flow rate to the orifice was varied from 0.3 
cm?/sec to 1.22 cm?/sec. giving a frequency range of 3 to 10 bubbles/ 
second. The ammonia-nitrogen feed mixture was kept at a constant comp- 
osition of about 19% ammonia by volume. The major variables for each 
run are listed in Table 12. 

The data are tabulated in Table 13 and the results are sum- 
marized in Table 14. Table 4 is a detailed data sheet for run number 
N2 and sample calculations for this run are presented in the appendix 
(A.4). 

The surface areas and volumes of the detached bubbles as 
obtained from the films are listed in Table 15 and Table 16 respect- 
ively. The values give a measure of the reproducibility of the ex- 
perimental technique. 

Figure 7 is a graphical representation of the variation of 
surface area of the bubble with time during the formation period. The 


variation is linear except in the neck formation zone. For some reason 
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Table 15: Surface Area of Detached Bubbles 


for NH3-N2-Water System 
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NO w SOE Volume of Detached Bubble, cm? 


Bubble 


Table 16: Volume of Detached Bubbles for NH3-N2-Water System 


<< on 


a tt Se 


on nue ¥ hy ows 
=— 2. 
{ 


a 


mo ,Slddvd Dads 


_ a) aie tine OO Ay pny aw om 


01.0 | a 
: 
01.0 | ee8o.0 
Aa 4 
Ladd .0 PYRO, 0 
ary 7 | ce Ri £ 
b. 
7 
1.0 | £ Bat 
cOl. BY ($#0.0 
Bt . HOCH |! 
: 
y\ bot - 


— 


Sf01,0 - 


/ 


j 
t | 
PCUL.\V : 


ie ee i i ee te ns 


ia 


10] 


ceo 


WAISAS UYSIVM—°N —*HN OS 
319gng ONIMOYNS Vv JO vaaVv “L:9l4S 


<= so3S * aWNIL 
820 v2°0 020 910 zl0 800 vOO 


WD 807 O0='AQ!l adAdIsI¥O 


V3aYvV J DVAYNS 


4 


MD 80>.0= a! 


85.0 “PSO 050 ato S10 


—— 20982 . Suit _ =F ae Ba 
| oa AS 
= 3Jj8GUG DMIWORD A WO ARAA VSI | ' 


MaTeve ASTAW-.4—,HM fOF = 


102 
the growth of the bubble is not well defined for run N2 indicating an 
unsteady growth of the bubble. The area of the rest bubble was added 
to the area of the detached bubble since at the instant of bubble 
detachment both the rest bubble and the detached bubble are contrib- 
uting to the mass transfer. 

The surface area of the bubble at any time during its growth 
could be approximately correlated to the frequency of bubble formation 


and time, from the data of Figure 7. The resulting relationship is, 


1.34 


A= (m/4) D.* + 0.507 0 too eee 22 ------------------------- (5.7) 
where, A = Surface area of the growing bubble at time 
"et" (cm?) 
Dy = Inside diameter of the orifice, (cm) 
n = Frequency of bubble formation, tnecu) 
t = Time at any instant of bubble formation, (sec.) 


The basis for arriving at this equation was the same as that 
for the carbon dioxide-water system but this equation is not strictly 
valid for the necking region of bubble formation. 

The volume of the growing bubble was also found to vary line- 
arly with time during the entire growth period. The results are shown 
in Figure 8. Once again the volume of the rest bubble was added to the 
volume of the detached bubble. 

All the data of Figure 8 could be approximately represented in 
the following form and the equation quantitatively shows the effect of 


frequency on the bubble volume at any instant of its growth. 
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V = V_ + 0.056 nit?? ,  ------------------------ (5.8) 
where, V = Volume of the growing bubble at time 't',(cm’) 
V_ = Volume of the bubble at time zero, (cm?) 


Oo 


Frequency of bubble formation, (seca) 


n 


tt 
It 


Time at any instant of bubble formation, (sec.) 


This equation is valid for the entire formation period of the 
bubble. 

The volume of the gas bubbles after detachment for the 
ammonia-nitrogen-water system was found to vary linearly with the fre- 
quency of bubble formation as shown in Figure 9. The range of frequen- 
cies studied was from 3 bubbles per second to 10 bubbles per second. 


Quantitatively, this could be expressed as, 


< 
i} 


10900 a+ 0200351 tiga een eee (5.9) 


(for 3 <n < 10) 


where, Ve = Volume of detached bubble, (cm?) 


Frequency of bubble formation, (seca 


n 


D Inside diameter of the orifice. (cm) 


fo) 
The gas phase mass transfer coefficient was calculated from 
the measured values of the overall transfer coefficient by correcting 


for the liquid phase resistance. The equation used is, 
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where, K. is the overall transfer coefficient, 


G 
gm-moles/hr-cm*-atm. 

k. is the gas phase transfer coefficient, 

gm-moles/hr-cm*-atm. 


k, is the liquid phase transfer coefficient, cm./hr. 


atm. 


m as the slope of the equilibrium curve, ———;—3 
gm—-moles/cm 


The value of 'm' at 1 atm. and 25°C is given in Perry's 


Handbook as 17.89 atm. /gm-moles/cm?* for the ammonia-water system. 


From the experiments done on the CO2-water system the value of 


the liquid phase mass transfer coefficient was measured to be about 50 


cm./hr. and it was found to be independent of the gas flow rate. This 


value of ke was assumed to apply to the NH3-N2-water system also since 


the diffusivities of ammonia and carbon dioxide in water are about the 


same (1.96 x 10F cem*/sec. for CO. in water and 2 x 10s cme eee, for 


NH3 in water). 


Thus equation (5.10) reduces to the following form for the 


NH3-N2-water system, 


or, 


teil cg L769 
Ko Ko 50 
re ss apes ecg te ne i wen a A i (5.11) 
oe: aa 
G 


(1-0.35275 Ko) 


The values of ke as calculated from equation (5.11) are plotted 


in Figure 10 and are tabulated in the appendix (Table 25). It can be 


observed that the gas phase mass transfer coefficient during bubble 


formation increases with the bubble frequency, or the gas flow rate. 
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CHAPTER 6 


DISCUSSION OF RESULTS 


6.1 Carbon Dioxide-Water System 


6.1.1 Volume of Detached Bubbles 

The volume of the detached bubbles for the carbon dioxide-water 
system, in the experimental runs where negligible amount of mass trans- 
fer was observed (Table 6), remained essentially constant at 0.094 c.c. 
for bubble frequencies of 3 and 5 bubbles/second showing that the 
bubble formation occurred in the constant volume region. Equation (2.4) 
is applicable in this formation regime as Reynold's number from definit- 
ion (Equation 2.1) is less than one in the range of flow rates studied, 
and predicts the value of the detached bubble volume to be 0.0945 c.c. 
when the appropriate values of variables are used. Thus there is a 
good agreement between the experimentally determined and the predicted 
values of the detached bubble volume. 

The final bubble volume was found to increase slightly with the 


bubble frequency. 


6.1.2 Area-Time Relationship For The Growing Bubble 
Calderbank and Patra [15] did an experimental study of carbon 
dioxide bubbles during their formation in water and found an empirical 
correlation for the area of the bubbles, which was given as, 
CA See 


= 1 -(4 11% ---------------------------------------- (2.9) 
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This equation was valid for bubble frequencies of 0.2 to 2.0 
bubbles/sec. and equivalent spherical bubble diameters of 0.5 to 1.0 
cm. In this equation A and AN are the area at any time t and the max- 
imum area of the bubble during its formation respectively. 

In Figures 11 and 12 a comparison has been made of the values 
of area predicted by equation (2.9) with the ones found from this ex- 
perimental investigation. The findings of Kalra [49] are also super- 
imposed on the graphs. In Figure 11 Kalra's results were obtained from 
a study of mass transfer during bubble formation using the carbon dio- 
xide-water system. The bubble frequency was 1 bubble/sec. and the 
water velocity in the column was constant at 1 cm./sec. His results in 
Figure 12 correspond to a frequency of 1/2 bubble/sec. and a water 
velocity of 1 cm./sec. The results of this investigation compare very 
well with Kalra's findings. However, the growth of the bubble was 
observed to be more steady in terms of its surface area in this study. 
In all the experimental runs the area of the bubble dropped after break- 
off from the orifice. This may be due to the stretching encountered 
by the bubble prior to its breaking which results in an elongated neck 
and hence a larger surface area. 

The equation proposed by Calderbank et.al. [15] does not agree 
with the findings of this study. They have not included the effect of 
the rest-bubble in their equation. Besides the equation predicts a 
very uniform increase of area with time till the bubble detaches. This 
is contrary to the results reported here and by Kalra [49]. In both of 
the latter studies a very sharp increase of area was noticed in the 


neck-formation zone of bubble growth. It can also be observed that 
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SURFACE AREA OF BUBBLE 
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Calderbank's equation overestimates the value of the area under the 
area-time curve. 

A new relationship for carbon-dioxide bubbles forming in water 
from a submerged horizontal orifice of 0.408 cm. is found to exist and 


is expressed in equation (5.4) as, 


2 


A = ma + 0.66 nt (for 'n' in the range of 0.5 to 1 bubble/sec.) 


4 


This equation is plotted in Figure 5 and it can be observed 
that the agreement with the experimental values is fairly good except 


in the neck formation region. 


6.1.3 Volume-Time Relationship For The Growing Bubble 

It was found that the volume of the bubble varied linearly with 
time in most of the formation region. The exception was in the initial 
period of the growth cycle where it grew a little slowly. This could 
probably be because of a sudden reduction in pressure inside the bubble 
after the bulk of it gets detached from the orifice. Since the water 
head is constant it would have a tendency to force the bubble inside 
the orifice. The linear variation of volume of the bubble subsequ- 
ently implies a steady gas flow rate to the orifice. Kalra [49] has 


also reported similar findings. 


6.1.4 Liquid Phase Mass Transfer Coefficient 
The value of the liquid phase mass transfer coefficient during 


bubble formation for the carbon dioxide-water system did not vary 
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appreciably with a superimposed co-current flow of the continuous 
phase. It held constant at about 45 cm./hr. at 1 bubble per second 
over a water velocity range of 0.3 cm./sec. to 2.5 cm./sec. A slight 
increase in the value of ke was observed at 0.5 bubble per second bubble 
frequency and the value remained constant at 49.5 cm./hr. over the same 
range of water velocity. Though the effect of water velocity was not 
isolated by Kalra [49] he found the values for a similar system to be 
160 cm./hr. at 1 bubble per second and 100 cm./hr. at 0.5 bubble/second. 
His study was done at a water velocity of 1 cm./sec. Kalra's values of 
the transfer coefficient during formation are higher than those observed 
by other workers for bubble rise period. 

Calderbank and Patra [15] have indicated that their data on a 
study of mass transfer during bubble formation using the CO2-water 


system agreed well with the equation, 


me 
k= 1-56 "Gee ---------- = -- $= = == $= = == = = $= $5 == = = === (27.27) 


where ky was defined as in equation (2.26). 

For the carbon dioxide-water system at 25°C the value of the 
GIfLUsivity is 1.96 x 10 cm?/sec. which gives the mass transfer 
coefficient at bubble formation as 14 cm./hr. for a bubble frequency of 
one bubble per second and 9.5 cm./hr. for half a bubble per second. 
These values are lower than the ones observed in this investigation. 

When Calderbank et. al.'s definition of kK is changed to 
equation (3.62) in order to make it consistent with the definition 


which is used to calculate ky in this study. the following equation is 


obtained, 
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which predicts the value of kK to be “21.6 °cms /hir? and’ 1456 
cm./hr. for bubble frequencies of 1 per second and 0.5 per second res- 
pectively. The values are still lower than the ones obtained exper- 
imentally by us. 

By modifying Calderbank and Patra's model to incorporate the 
area-time relationship of the form defined by equation (3.57), it was 
derived that, 


Dy A, ret) 670k" t, 
k= 2 "Geol, 0.5 kre, re ATE eons re ae ea eas (3.60) 


From the experimental results k' could be expressed as a 


function of the bubble frequency as in equation (5.2), 


which on substitution in equation (3.60) gives, 


A + 0.442 
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Further, the surface area of the rest bubble was found to remain 
essentially constant and could be approximated by equation (5.3). In 
this study the internal diameter of the orifice used was 0.408 cm. and 


thus equation (6.2) reduces to, 
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and it predicts the values of the mass transfer coefficient during 
formation to be 22.2 and 15.1 cm./hr. for bubble frequencies of one and 
one half per second respectively. Although these values are low as 
compared to the experimentally observed values of 45 and 49.5 cm./hr., 
they agree very well with the predictions of Calderbank's modified 
equation (3.63). 

Beek and Kramer's [3] solution was shown in equation (3.45) for 


a spherical bubble growth to be, 


The values of ke calculated from this equation for the carbon 
dioxide-water system are 17.9 cm./hr. at one bubble per second, and 
12.2 cm./hr. at half bubble per second. The values are again lower 
than the experimental measurements but are of the same order of magnit- 
ude as those predicted from equations (3.63) and *( 6:43). 

Popovich et. al. [ 67] have proposed an equation of the form 
(2.23) for mass transfer during drop formation and found tke following 
equation to be valid from their experimental findings on the sodium 


iodide drops formed in isobutyl alcohol. 
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Defining, 


where, A for spherical drop growth is given from equation (2.24) 


as, 
. 2 -2/3 .2/3 
A=T d_ Ce i he Spach aaa age a ig wake od men ee ae gl ale eel J cle te ae (6.6) 
From equations (6.4), (6.5) and (6.6) it can be shown that, 
D 
4 Bere oe eS a es ee 
ioe= 2.2 "GED -------------------- (6.7) 


This equation, once again, predicts the value of kK in the same 
range as obtained from equations (3.63), (3.45) and (6.3). 

The mass transfer during the rise period from a circulating 
bubble can be predicted from equations (2.20) and (2.21) when the 
Reynold's number is greater than 400. The detached volume, V, for this 
Study was found to remain constant at about 0.077 c.c. from Table ll. 
Hence, from equation (2.21) the terminal velocity of carbon dioxide 
bubbles in water is calculated as 7.6 cm./sec. which is very low when 
compared to the corresponding experimental values listed in Table l. 
The volume of 0.077 c.c. gives an equivalent bubble diameter of 0.528 
cm. for which the velocity of bubble rise is about 20 cm./sec. 

With this value of velocity the Reynold's number from equation 


(2.10) is 1055 and hence equation (2.20) can be used. It predicts Ke 
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kz 
during rise to be 110 cm./hr. which is about the same as that observed 
experimentally by various investigators listed in Table l. 
A comparison of the above value of transfer coefficient with 
the values obtained for the mass transfer coefficient during bubble 


formation in this study shows that at a bubble frequency of 1 bubble/ 


second, 
k 
Mma tO) BS cs OU La eee eal ee ey OU (6.8) 
eT (rise) pe 
and at half bubble per second, 
bt (f tion) 49.5 
pate ele Ue ormation i Nl ic = O F 45 Pp ee hn era ene ae eee a gece a a ee ne (6 . 9) 


eT, (rise) 110 


In article 3.2.1 by using a surface stretching model for the transfer 


coefficient during formation and Higbie's penetration theory for ko 


during rise it is shown that, 


ke (formation) d, 
kau(fise) = Ne an race ae epaaIais Gas ae ee Oe (3.50) 
where, 2 = U ee ey EES OR TT a (2.29) 


As before, d. in this study was approximately 0.528 cm. and from 


a 


Table 1 the corresponding value for U is 20 cm./sec. Thus equation 


(3.50) simplifies to, 


k -(0.5 
ple (hotmat torn yee Cali ante, Dubble forme liens 2 (6.10) 
My (rise) 


bevreedo Jedi 26 omse-ostls sucda ef askew. 
es 


.! aldsT nt batelti etosugises } 
“ : L 
djtw jusrstiiecs astansyy 36 ‘ebley avails sila do" 


olddud_saiisk jnsiotise7 1otenas3 sasig-ssld. x02 ‘eeaal 
pe a = zi _ 
\siddud I io vageupsy? siddbd & 36 Jedd awoile Nivap 4 et a) 


~, 


» 
a4 
nN 
a | 
nes ae eens meee MTEL 
(v +} ——e- — oS tae eee “7s .f) 


tsiengid sii ato? lebom grirptseyvie onaiwe 


A. 10% yrosds nOltsxtsnsg Sr -SidyiH itis no sama nto 

iv ore > 

, Jens wo at 
saya 
aginiy 


I y , 


fe 

(O25) oe ae a ee eee a 

(@¢,S laced ieetineiatiatiieeieel ll cetieeieatenl online lanlententeelalentantentedienhiientiendamamemmninetemamiall -U. 7 
) 74 7: 


> Arne 


‘ an, oe 
met? bone .mo 8S¢.0 yistiatxorgqe. aby siatva: ales nk y bt 


noljsups aval sda. sir vu 103 ouley 3 


118 
which gives the ratio as 0.162 at 1 bubble per second and 0.115 at 1/2 
bubble per second. These values are much lower than those observed 
experimentally (equations 6.8 and 6.9). 

The model described in article 3.2.2 predicted the ratio of mass 
transfer coefficients during bubble formation and rise as equation 
(3.64). This equation on comparison with equations (2.29) and (5.4) 
simplifies to the following equation when the diameter of the orifice 


is taken as 0.408 cm. 


k (formation) d 
kK, ise) ml. 24 vas) a lon 9 eal ear = nee rnd lel lan nr (Grit) 


Thus the ratio predicted by equation (6.11) is 24% higher than 
those by equation (3.50) or equation (6.10). The values of kK (form- 
ation) /k, (rise) are 0.2 and 0.143 for phe gubacias of 1 and 1/2 bubble 
per second respectively. Once again, these values are smaller than the 
ones in equation (6.8) and (6.9). 

Hence the mass transfer coefficient during the bubble formation 


is lower than its value during the subsequent bubble rise period. 
6.2 Ammonia-Nitrogen-Water System 


6.2.1 Volume of Detached Bubbles 

The volume of the gas bubbles after detachment for the 
ammonia-nitrogen-water system increased linearly with the frequency as 
in equation (5.9) showing thereby that the experiments were done in the 


slowly increasing volume region of bubble formation. 
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6.2.2 Area-Time Relationship for the Growing Bubble 

The area of the bubble increased linearly with time but varied 
as iva es power of the frequency of bubble formation. The relation- 
ship is expressed in equation (5.7) and is different from the one ob- 
tained for the carbon dioxide-water system. However, the necking region 
is not well represented by equation (5.7). The variation in the surface 
area of the rest-bubble was not much (0.11 to 0.13 cm? ) for the 0.408 
cm. internal diameter orifice used and in equation (5.7) it is approx- 


imated to 0.13 cm’. 


6.2.3 Volume-Time Relationship for the Growing Bubble 

The volume-time relationship was found to be much better for the 
ammonia-nitrogen-water system than for the carbon dioxide-water system. 
In the latter system a deviation from the linearity existed in the in- 
itial growth period of the bubble. On the contrary, no such effect was 
observed for the former system. The volume of the bubble varied as 
frogetl power of the bubble frequency whereas for the CO2-water system 
it varied as a linear function of the bubble frequency. The volume of 
the rest-bubble in equation (5.8) remained essentially constant and had 
a value of 0.0696 cm*® for the 0.408 cm. internal diameter orifice used. 

The linear variation of bubble volume with time gives an in- 


dication of the steadiness of the gas flow rate to the orifice. 


6.2.4 Gas Phase Mass Transfer Coefficient 
The individual gas phase mass transfer coefficient during 


bubble formation was found to increase with the gas flow rate to the 
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orifice and the variation is depicted in Figure 10. The value of the 
transfer coefficient was in the range of 0.25 to 0.32 gm-moles/hr-cm?- 
atm. over a bubble frequency range of 3 to 10 bubbles/second. The 
orifice inside diameter was 0.408 cm. and the composition of the 
ammonia-nitrogen feed was about 19% ammonia by volume. The mass trans- 
fer during the formation period was such as to decrease the bubble vol- 
ume by 13 to 17 percent. 

Sharma and Mashelkar [73] have determined the values of gas- 
side mass transfer coefficients during bubble rise by absorbing sulphur 
dioxide bubbles diluted with air in aqueous solutions of sodium hydr- 


oxide. Their experimental value of k, was 0.203 gm-moles/hr-cm*-atm. 


G 
The corresponding value for ammonia is estimated to be 0.281 
gm-moles/hr-cm*-atm. This correction is based on the assumption that 
ko is directly proportional to the square root of the diffusivity of 
solute in the gas phase. The diffusivities for SO2 in air and NH3 in 
Np are 0.133 cm*/sec. and 0.254 cm*/sec. respectively. Thus the values 
observed experimentally by us for ke during formation are of the same 
order of magnitude as those obtained during bubble rise. 


Geddes' equation (2.42) predicts the value of k, during rise 


G 


from a stagnant sphere and for a:gas contact time of t, reduces to, 
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and for the ammonia-nitrogen feed the value of De ac. 25: Cand’ 1l-atm.. is 


0.254 cm*/sec. which gives, 
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Serato) = (9-165 R + eo cm./sec. amen nnn aan ae ee (6.13) 
t R 4 
ft 
vee KG (rise) = 0.147 Chai Lat a (in Sancta h --- (6.14) 
t. R hr-cm -atm. 


This equation on substitution for the value of bubble radius, 
R, from the experimentally determined vclumes of the detached bubbles 
shows Ko to be 0.47 gm-moles/hr-cm*-atm. at 3 bubbles/second and 0.48 
gm-moles/hr-cm*-atm. at 10 bubbles/second. These values of ke are 
higher than those observed experimentally (0.25 and 0.32 respectively) 
during bubble formation. 

Calderbank's equation (2.38) also predicts about the same 
values of ko as equation (6.14). But this is expected since the only 
difference in the two equations is the term 0.165 R/t,, which is small 
in comparison to the term De Te L3Ry 

Vivian and Behrmann [83] have studied the absorption of amm- 
onia into distilled water in a short wetted-wall column for various 
compositions of an ammonia-nitrogen gas mixture. They found the value 
of the gas phase mass anstes coefficient to be about 0.2 gm-mole/ 
hr-cm*-atm. for an ammonia concentration of 19% by volume. This value 
is of the same order of magnitude as those found experimentally in this 
investigation for the bubble formation period. 

The gas phase mass transfer coefficient during bubble formation 
can be predicted by Thorogood's equation (2.48). This equation for the 


ammonia-nitrogen feed at 25°C and 1 atm. simplifies to, 
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where n is the frequency of bubble formation, sec 

A comparison of the values obtained from equation (6.15) with 
the experimentally measured values shows the latter values to be 
higher than the predicted values. At 10 bubbles/second the experimental 
value of ke is 0.32 gm-moles/hr-cm?-atm. whereas the calculated value 


from equation (6.15) is 0.26 gm-moles/hr-cm*-atm. and at 3 bubbles per 


second these two values are 0.25 and 0.15 respectively. 
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CHAPTER 7 


CONCLUSION 
7.1 General 

The experimental method used in this investigation is suitable 
for studying the mass transfer phenomenon during bubble formation pro- 
vided the amount of mass transfer occurring is such that the bubble vol- 
ume decreases by 5 percent or more. This limitation is due to the ex- 


perimental uncertainty involved which was + 22. 


7.2 Carbon Dioxide-Water System 
The value of the liquid phase mass transfer coefficient during 
bubble formation for the carbon dioxide-water system is lower than its 
value during bubble rise and the ratio of the two coefficients is about 
0.43. The theory predicts a ratio of about 0.2 which is even lower 
than that observed experimentally. 
ky during bubble formation is independent of the liquid velocity 
and the bubble frequency over the range studied. The value of ky 
measured is about 47 cm./hr. over a water velocity range of 0.3 to 
2.5 cm./second and a bubble frequency range of 0.5 to 1.0 bubble/second. 
The experimental value of the transfer coefficient can be used 


along with the equation, 


to estimate the amount of mass transfer occurring during the formation of 
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a bubble under the liquid phase controlling conditions. The area of 


the bubble is given by, 


For example, by using this method the amount of mass transfer 
occurring at 3 bubbles/second is estimated to be 0.0052 cc/sec. which 
is in fair agreement with the experimental values. 

In additon, a rough approximation of the transfer coefficient 


is given by the modified penetration theory as, 


D 0.786 DG + 0.442 


Meee cee ogee? +0.33) 0 ye (6.2) 


7.3 Ammonia-Nitrogen-Water System 


The individual gas phase mass transfer coefficient during bubble 
formation was found to increase with the gas flow rate to the orifice. 
The value of the transfer coefficient was in the range of 0.25 to 0.32 
gm-moles/hr-cm?-atm over a bubble frequency range of 3 to 10 bubbles/ 
second. The water velocity in the bubble column was kept constant at 
1 cm./sec. and the composition of the armonia-nitrogen feed was maint- 
ained at about 19 percent ammonia by volume. 

The gas-side transfer coefficient during bubble formation is of 
the same order of magnitude as that during bubble rise. 

When the gas phase resistance is controlling the mass transfer 
occurring for a growing bubble can be estimated by using the exper- 


imental value of the coefficient into the equation, 
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Where (P.-P;) on is the logarithmic mean driving force during 
the formation period, and A is the area of the bubble at any time and 


is, 
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NOMENCLATURE 
Area, cm?, 
Interfacial area per unit volume of the dispersed phase. 
Point source (= w/4m), cm?/sec. 
Surface area of the bubble, cm”, 
Maximum area of the bubble during formation, cm’, 
Surface area of the rest bubble at zero time, em’, 
Area at any time t. 
Proportionality constant in equation (2.50). 
Acoustic velocity in the gas. 
Concentration, gm-moles/cm?. 
Concentration of solute in bubble after time te. 
Solute concentration at time t. 
Bubble diameter. 
Equivalent spherical diameter of bubble. 
Final equivalent bubble diameter, cm. 
Column diameter. 
Effective dispersed phase diffusivity. 
Diffusivity in the gas phase, cm7/sec. 
Diéfuetvity in the liquid phase, cm*/sec. 
Orifice internal diameter. 


Fractional approach to equilibrium ey, 
Height of bubble top from the orifice are at any time t 
Total amount of gas absorbed, gm-moles. 

A constant. 


A constant. 
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Slope of area-time curve for growing bubble, em?/sec. 
Slope of volume-time curve for growing bubble, cm? /sec. 
Gas phase mass transfer coefficient. 
Liquid phase mass transfer coefficient. 
Overall gas phase mass transfer coefficient, em-moles/hr.-cm?-atm. 
Overall gas phase mass transfer coefficient, cm./hr. 
Vertical distance between rising bubble centres. 
Mass transfer rate, gm-moles/sec. 
Mass transfer during bubble formation, gm-moles/sec. 
Mass transfer flux during bubble formation, gm-moles/cm?-sec. 
Frequency of bubble formation, geatln 
Rate of diffusion, gm-moles/cm?-sec. 
Mass transfer per unit area, gm-moles/cm?. 
Total mass transfer during bubble formation, gm-moles. 
Capacitance number. 
Partial pressure in the bulk of the gas. 
Polar radius. 
Orifice outside radius, cm. 
Equivalent spherical radius. 
Bubble equivalent radius, cm. 
Final equivalent bubble radius, cm. 
Surface area of the growing bubble, cm*, 
Surface area of the rest-bubble, cm’, 
Time. 
Exposure time. 


Time of formation of bubble, sec. 
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Tangential velocity component. 

Velocity of rise of bubble, cm./sec. 
Velocity of gas at the orifice. 

Bulk liquid velocity, cm./sec. 

Radial velocity component. 

Volume of bubble, cm?, 

Chamber volume. 

Volume of the bubble at zero time, cm?, 
Volume of detached bubble, em? , 


Volume at any time t. 


Constant characteristic of the system. 


Distance into the bulk phase normal to the interface. 


Distance measured in the direction of diffusion. 


Greek Symbols 


B 
6 


p(t) 


Constant in equation (2.23). 

Film thickness. 

Polar angle. 

Viscosity 

Kinematic viscosity, centipoise. 

Density. 

Surface tension, dyne/cm. 

Rate of mass transfer at time t, gm-moles/sec. 


Gas feed rate to the orifice, em?/sec. 
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Dimensionless Numbers 


Pe = Peclet number, d,U/D, . 

Re = Reynold number, dU of Hee 
Resi Orifice Reynold's number, Py Uy Difu ot 
Sc = Schmidt number, u c/o De 
Sh = Sherwood number, k d,/D, . 
Subscripts 

A = Component A. 

c = Continuous phase. 

d = Dispersed phase. 

e = Equilibrium values. 

g = Gas phase. 

a = Conditions at interphase. 
L = Liquid phase. 


O = Initial condition. 
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APPENDIX 


-L Calibration of Liquid Flow Rotameters 


Two rotameters R13M-25-3 and R-10M-25-3 were calibrated to 
measure liquid velocities ranging from 0 to 2.75 cm./sec. in a column 
of 20.2 cms. internal diameter. 

The calibration data is given in Table 17 and the calibration 


curves are shown in Figure 13. 


.2 Calibration of Gas Flow Rotameter 
Matheson flowmeter No. 610 was calibrated for air and ammonia 
separately at 1 atmosphere and 70°F. The calibration curve is shown 


in Figure 14. 


A.3 Calculations for Determining ko from CO2-Water Data 


As an example the data in Table 3 was used for run M8 and 
the calculations would be done stepwise. 

1. The time period data was analysed by a computer program to 
calculate the sample size. This sample size was used to det- 
ermine the number of bubbles to be analysed for volume and 
surface area on the digitizer. 


The statistical formula used was, 
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Average Liquid Velocity 
Rotameter Reading In Column, cm./sec. 


R1OM-25-3: 
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REOM<-253 : 
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Table 17: Liquid Flow Rotameters Calibration Data 
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where, N is the sample size 
n is the number of data points 
x is the arithmetic mean time period 
Xs is the fie time period 
It gives a 99% confidence limit. 
e.g. in this case N = 13. 
It was decided to analyse 10 bubbles. 
The average recorded flow rate of the gas was converted to the 


flow rate at orifice conditions by using the ideal gas law. 


15.52 (PSIA) 5 Meee (°K) 
14.96 (PSIA) ~ 295.4 (°K) 


Flow at orifice x 0.0822 c.c./sec. 


0.0860 c.c./sec. 


Ten detached bubbles were analyzed on the digitizer and the 
volumes and surface areas calculated numerically by using 

a computer program written by Kalra [49] and the mean bubble 
volume was found. 

The bubble having a volume closest to the mean value 
was analyzed at various points during its formation and the 
calculated area from the computer was plotted as a function 
of time as in Figure 5. 


The area under this curve was found to be, 


if f A dt = 0.446 cm?-sec. 
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The CO2 concentration at the interface was calculated from 
the solubility data. 

Temperature of water = 25°C. 

Assume Henry's law to hold for CO2-water. 
3 


Henry's constant, H = 1.64 x 10. 


Also, Cc, = Pa 


moles CO2 


where, C, = CO2 concentration at the interface, 
i mole soln. 


Pa 


partial pressure of CO2, atm. 
The gas COz2 enters the orifice saturated with water 


vapour and hence the partial pressure of CO2 is equal to the 


pressure at the orifice minus vapour pressure of water at 25°C. 


Pa = (14.96 - 0.46) psia 


= 14.5 psia 
= 0.987 atm. 
C ba 0.987 gm-moles CO2 


1.64 x Wits Papago cece solution 


0.334 x 107 gm-moles CO2 


c.c. solution 


CO2 concentration in the bulk of the liquid was taken as zero 
since 12.6 c.c. of acid was required for both the blank and 
water sample. 

Mass transferred was calculated from the mean bubble volume, 


time period and gas flow rate at the orifice. 
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Step 4. 


Mean bubble volume = 0.077 c.c. 
Mean time period = 0.9604 sec. 


Gas flow @ orifice from Step (2) = 0.086 c.c./sec. 


0.077 


(0-086° = 5 060d) 


cee. (sec... 


B 
i} 


8.94 x 104 (gm-moles C0O2/hr.) 


Average liquid phase mass transfer coefficient, 


k, SS pe = REE cme ir. 


£ 
oi A dt) C, 


e13.94° x 107" x 0.9604 


a em./hr. 
0.446 x 0.334 x10 


57.6 cm./hr. 


Calculations for Determining ko From NH3-N2-Water Data 


Sample data sheet shown in Table 4 for run N2 would be 


used. 
Steps 1, 2, 3 were done in the same manner as for CO2-water. 


Flow at orifice = 0.3026 c.c./sec. 


M A dt = 0.189 cm?-sec. 


Taking a mass balance for ammonia gives, 
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(NH3 flow in + NH3 in rest bubble) 
= (NH3 in detached bubble) + (NH3 in rest bubble) 


+ mass transfer 


Symbolically, 
Coy vit nem toe tye y 
: at rae t be Jee 
where, G = Gas flow rate @ orifice, c.c./sec. 
Fie Inlet ammonia mole fraction in gas 


Ye = Final ammonia mole fraction in gas 


v, = Volume of detached bubble, cm? 
m, = Mass transfer, c.c./sec. 
te = Formation time, sec. 


weve Oo Ys = N,)/(v_/t,) 


Assuming bulk liquid concentration to remain sub- 
stantially zero. 


Driving force for mass-transfer during bubble formation, 
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m  An(y,/y,) 
(0.3026 - ee) 
_ 0.3026 x 0.184 - ie 
here, Ye = 
(0.0805/0.3231) 
= 0.00932 
y_ = 0.0586 
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Step 5. Average overall 


“alae 

Or; Koc = 

Also Ko = 
where, Ko A 
R is 

i is 


Hence at 25°C., 


A.5 Bubble Area and Volume vs. Time Data for the CO2-Water System 


gas phase mass transfer coefficient, 


meas 


8 


f 
Cf 


A dt) NA 


0.0534 \(c.¢c./sec.) x 0.3231 €sec<) 


0.189 (cm*-sec) x 0.0586 


1.557 em./sec. 


= 5610 cm./hr. 


Koc /RT 


in (gm-moles)/(hr) (cm?) (atm. ) 


82.06 (atm. ) (cm?) / (°K) (gm-mole) 


0 
1ty F ek 


Ko = 0.229 (gm-moles)/ (hr) (cm?) (atm) 


The data are presented in Tables 18 to 24. 
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Time During 
Bubble Surface Area, cm? Volume, cm* 
Formation, Sec. 


Table 18: Bubble Area and Volume vs. Time for Run M6 
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Bubble Surface Area, cm? Volume, cm? 
Formation, Sec. 


Table 19: Bubble Area and Volume vs. Time for Run M7 
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Time During 
Bubble Surface ‘ Volume, cm 
Formation, Sec. 
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Table 20: Bubble Area and Volume vs. Time for Run M8 
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Table 21: Bubble Area and Volume vs. Time for Run M9 
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Time During 5 ' 
Bubble Surface Area, cm Volume, cm 
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Table 22: Bubble Area and Volume vs. Time for Run M10 
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Surface Area, cm Volume, cm 


Table 23: Bubble Area and Volume vs. Time for Run M1l 
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Table 24: Bubble Area and Volume vs. Time for Run M12 
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A.6 Mass Transfer Coefficients for the NH3-N2-Water System 


xe Kg 
gm-moles/hr-cm?-atm gm-moles/hr-cm?-atm 


Table 25: Mass Transfer Coefficients for the NH3-N2-Water System 
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